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CHATTER  III 

Problems  of  Estimation 

It  ie  frequently  important  that  one  sake  estimates  of  mean  life, 
rates  of  failure,  probability  of  survival  for  a  given  tine,  eta,  on  tike 
basis  of  data  arising  fras  life  tests.  The  data  say  be  generated  in 
■soy  ways;  e.g. ,  they  nay  arise  from  truncated,  censored,  sequential, 
replacement,  noa-replaceacat,  interrupted,  or  combined  experiments; 
we  nay  or  nay  uot  know  the  exact  tines  to  failure.  We  shall  try  in 
what  follows  to  give  rules  and  procedures  which  enable  us  to  give  point 
and  interval  estimates  which  axe  in  some  sense  optimum. 


vt* 


I 

i 

i 
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Section  1. 


Estimation  in  the  Censored  One  Semple  Case.  (Humber  of  f alltires  is 
fixed.  Items  which  fail  r-#y  or  may  not  he  replaced). 

Basic  Considerations,  fra  lot  and  Interval  Estimates  for  6  . 

Let  ua  make  the  following  assumptions: 

(i)  n  1  teas  are  drawn  at  randan  from  a  density  function  of  the 

form  f(x;0)  »  ^  e“x^,  x  >  0,  8  >  0; 

(11)  the  n  items  era  placed  on  life  test  at  time  sero  and 

failure  tines  become  available  In  order.  That  is  to  say, 

X,  _  <  x.  <  . . .  <  x  <  . . .  <  x  ,  where  by  x.  ^  ie  meant  the 
l,n  —  <!,d  -  —  r,n  —  —  n,n  x,n 

time  when  the  itb  failure  occurs,  (measured  from  the  beginning  of  the 
life  test). 


v 

4> 


I 


) 


i 


(Hi)  the  experiment  is  discontinued  as  soon  as  x  has  he  corns 

r 

available  (i.e.,  after  the  first  r  observations  are  made). 

He  wish  under  assuaptionu  (i),  (ii),  and  (ill)  to  find  a  "good" 
estimate  of  9  and  to  give  the  distribution  of  this  estimate  in  both 
the  non- replacement  cane  (where  failed  Items  are  not  replaced)  and  in 
the  replacement  case  (where  failed  items  are  replaced  immediately  by 
new  items).  This  is  given  by  the  following  theorem: 

Theorem:  Under  (i),  (ii),  end  (ill)  an  estimate  based  on  the 
first  r  ordered  observations  which  is  "best"  in  the  sense  that  it  is 
maximum  likelihood, unbiased,  minimum  variance,  efficient,  and  sufficient 
is  given  by 


(1) 


A 

6 


r,a 


T„/r 


where  T^  is  the  total  life  of  items  on  test  observed  up  to  the  time 
of  the  r01  failure.  In  the  non- replacement  case: 


(2)  Tr  M  wtk  +  (n~l)(x2-xI )  +  ...  +  (n-i+l)(xi-x1_1)  +  ....  +  (n-r+l)(xr-xr  x) 


■  XI  xi  +  (“-r)xr  > 


i«l 


and  so  the  "best"  estimate  ( i)  becomes 


(3) 


e 


,n  ’  |g>l  *  !""r)jS- 


/'  • 


■r 
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In  the  replaceaent  cnee: 

(4)  *r  -  nx^  +  nfXg-J^)  +  ...  +  n(xr-xrl)  *  n*r 

and  so  the  "best"  estimate  (1)  becones 

(5)  ®r,a  -  oxr/r  . 

A 

The  probability  density  function  of  0  in  either  the  replace^ 

r,n 

meat  or  non-replacenent  case  is  given  by 

(6)  fy(y)  ■  (-rryt  (r/e)ryr'’1e‘ry/e,  y  >  o 

a  0  ,  elsewhere. 

Ibe  proof  of  this  theorea  is  given  is  Appendix  3B. 

Pro*  (6)  it  follows  at  once  that  W  «  2r  6^  a/e  *  2Tr/6  is 
distributed  a a  ^(2r).  Consequently  if  the  constant  X*(2k)  is 
defined  as  Pr(X?(2k)  >  X^(21i)  .«?  y,  then  a  lOO(l-a)  percent  two- 
Glded  confidence  Interval  for  6  is  given  by 
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f 

V 

& 

«. 

I 

ft' 


P 


f 


2T 


equivalently 


.xa(^) 


v’“)t 


will  cover  the  true  hut  unknown  value  of  9,  lOO(l-a) 


percent  of  the  tine. 

Let  c^(r,a)  m  2r/*^(2r)  ,  then  a  lOO(l-o)  percent  one-sided  confidence 
interval  for  6  can  be  written  ae  [c,(r,a)e  ,  oo )  .  In  Table  2  we  give 

J  *  )“ 

the  values  of  c^(r,a)  for  a  m  .01,  .05,  .10,  .20,  .25,  and  .50  and 
r  <•  1(1)  20(5)  30(10)  50(25)  100  . 

For  large  r  (say  >50)  >^(2r)  ia  approximately  normally  ciutributed 
with  mean  2r  and  variance  Vr  .  Consequently,  the  two-sided  loo(l-a) 
percent  confidence  interval  becomes  (for  large  r) 


(10) 


e 

~ca  ' 
1  +  — 


A 

e 


1  - 


where 


ca  «  2.576 
®  1.960 
m  1.645 
«  1.262 
*i  .67^ 


a 


if 


a 


•  .01 
■  „05 
*  .10 
*  .20 
*  .50 


ft 


In  the  one-sided  case  ths  lOO(l-a)  percent  confidence  intern  al  becomes 


(11) 


§ 


1  + 


^  * 


00  1 


;,v 

li. 


1 


*3'6* 


vhert 


dQ  *  2.326 

•  1.645 

-  1.282 

*  .674 


o  «  .01 

•  .05 
«  .10 
-  .25 

*  .50 


Estimation  of  Other  Quantities: 

(a)  In  manor  practical  problems  one  does  not  wish  to  find  point  or 

Interval  estimates  for  the  mean  life  6,  but  rather  for  a  quantile  x^, 

where  x  Is  that  life  such  that 
P 


Pr(X  >  *p)  •  P  • 


For  the  exponential  p.d.f.  this  swans  that 


-x/a 


«  p  or  X^  m  e  log  X 
p  p 


It  is  therefore  clear  that  the  maximum,  likelihood  estimate  of  x 

P 

Is  given  by  $  log  ~  ■  Furthermore,  tvo-slded  and  one-sided  lOO(l-a) 

percent  confidence  intervals  for  x  are: 

P 


A  1  <*  1 

2r0  log  ±  2r6  _  log  ~ 
— L£~-l  .  ,.?J? 


or  equivalently 


PB#**'**  «**>-.».  .  wwwwt'**  **#»  yr  wy  jnaary  yiBj^^WWPpagap) 
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s 


and 


(Hi) 


2Tr  log  |  gfr  log  |  | 

?^r  '  a<*>  / 

2  a'  ■* 


I  2r$  log  |  \  /  2Tr  log  §  \ 

- *  ,  03  or  |  -5= - *  >  00 

^  (to)  I  \  >§to)  I 


respectively . 

La  Table  3  we  give  values  of  log  -  for  various  useful  values  of  p  - 
Two-sided  and  one-sided  confidence  intervals  for  can  be  found  by  using 

Tables  1,  2,  and  3  and  substituting  appropriately  in  equations  (14)  and  (15). 
Remark  1:  Formula  (15)  can  be  interpreted  as  follows.  On  the  basis  of 

y4 

the  estimate  8  we  can  be  100(1  -or)  percent  confident  of  the  assertion 

r  jn 

that  the  probability  of  surviving 


2r6 


jog  ~  time  units 


£ 


is  >  p 


This  is  a  tolerance  Interval  statement  in  the  sense  that  if  we 


observe  0^  o  for  a  sample  we  can  be  100(1-0 )  percent  confident  of  the 
correctness  of  the  assertion  ramt  the  fraction  of  lteas  in  the  population 
surviving  x  or  mot'e  time  units  is  >  p  . 


Remark 


It  should  also  be  noted  that  if  we  observe  6  then 

-  ,n 


formulae  (14)  and  (15)  give  one  and  two-sided  100(1 -a)  percent  confidence 
bands  for  the  entire  distribution. 
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(b)  Frequently  v«<  vish  to  aake  confidence  stateaents  about  the 
proportion  of  Itens  surviving  some  preae signed  tie*  t#,  on  the  heals 
of  the  first  r  failure  tines .  Since  the  probability  of  surviving 
for  a  fined  tine  t*  is  given  by 


(16) 


Pt*  -  Pr(X  >  t*)  -  e't#/a 


it  is  clear  that  the  aaaclsua  likelihood  estieate  of  p^*  is  given  by 


(17) 


-t*/e 

.*  •  e  •  r,i 


Tram  (7)  it  follows  insediately  that  a  lOO(l-a)  percent  two-sided 
confidence  interval  for  pfc*  is  given  by 


-^(2r)t*/2rer>r 


1-  £ 


(18) 


or  equivalently 


\  e 

One-sided  100(1 -or)  percent  confidence  intervals  for  e  ®  are  particularly 
important.  It  is  an  iaasdiate  consequence  of  (9)  that  this  confidence  inter¬ 
val  is  given  by 


^(2r)t*/2Tr  -X2  a(2r)t*/2T. 

5  .  l'S 

9  e 


\ 
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I  */2re_  „  j 

(19)  U  r**  ,1 ) 


or  equivalently. 


•  \  *  » 


Die  question  nay  be  acised:  how  large  should  the  observed  0  (or 

r^u 

equivalently  T  )  be  in  order  that  ve  be  lOG(l-ct)  percent  confident 
that 


P^,®  »  e 


-**/«>  V  ? 


From  (19)  this  icplies  that 

(20) 


-<(SrK V2r»r  n 

e  ’  >y 


This  is  equivalent  to 

eT#n  >  ^(2r)t*/2r  log  y  or  Tp  >  ^(2r)tc/2  log  i  » 

The  naaning  of  the  Inequality  is  as  follovB; 

If  the  total  life  observed  In  getting  r  failures  exceeds 

y.2(2r)t*/£  log  4  ,  then  ve  can  be  100(l-a)  percent  confident  that  the 
G#  Jf 

probability  of  surviving  the  time  t*  is  >  y  .  These  values  are 
readily  confuted  from  Tables  2  and  3. 

Numerical  Examples 

Note:  It  is  assumed  throughout  that  the  underlying  distribution  of 
life  is  exponential. 


Vi 


ff.litUK-,  • 
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Example  1:  20  electron  tubes  are  placed  on  teat.  A  tube  which  fails 
is  replaced  at  once  Ij  e  new  tube.  The  fifth  failure  ia  observed  to  occur 
**07  hours  after  the  start  of  the  life  test. 

(a)  Estimate  the  mean  life  0  and  give  one  and  two-sided  9% 
confidence  intervals  for  9  . 

(b)  Estimate  x  _  ,  where  x  Q  is  such  that 

•y  *y 

! 

Pr(X  >  x>9)  -  .9 


Give  one  ««i  two-sided  9 5#  confidence  intervals  for  x  ^  . 

(c)  Make  &  one  and  two-Bided  95#  confidence  statement  for  the 
probability  of  surviving  100  hours. 

Solution: 

(a)  We  are  dealing  with  a  replacement  situation  with  n  *  20  ,  r  *  5  , 
x*.  m  U07  .  The  total  life  observed  is  given  by  T*.  ■»  20x,.  »  20(407)  *  8l4o  . 
Thus  it  foil  ovs  from  (3)  U..*i 


6  m  Tj/5  "  1628  . 


To  find  a  two-sided  95#  confidence  interval  we  use  (7)  with  “  20.483 


and  X^975(lo)  «  3-247 


This  gives  the  two-sided  interval  (T95>50l4)  .  To 


find  a  one-sided  95#  confidence  interval  we  use  (9)  with  X2  (10)  *  18.307  . 

»\jy 

This  gives  the  one-sided  interval  (809,co )  .  The  values  can  also  be  obtained 
directly  from  Tables  1  and  2. 


3-« 


;  I 


i  ■ 

(b)  Tbs  solution  is  found  by  Multiplying  through  by  log  -  «  log  ip  -  .103V. 
Thu*  vfe  get  -  (l628)(  .105^)  -  172  . 

A  y%  two-sided  confidence  interval  is  given  by  (83.8,528)  end  a  95)6 
one-sided  confidence  interval  is  given  by  (93*7>®  )  « 

(c)  The  maximal  likelihood  estimate  of  pt#  ,  the  probability  of  surviving 
t*  •  100  hours  le  given  by  pt*  «■  e~(lOO)/l628  _  e-.06l4  ^  9Jto4  similarly 

a  two-sided  95f  confidence  interval  for  pt#  is  given  by 

(,-100/TW  .-10O/5OH,  .  (.-1258,  .-0199,  .  (.881, ,.,802)  and  a  one-sided 

9 yf>  confidence  interval  for  pt*  is  given  by  (e_100^®^,l)  •*  (e~’'l42^,l)  » 
(.9936,1)  . 

Example  2:  20  electron  tubes  are  placed  on  test.  Tubes  which  fail 

are  not  replaced.  The  first  five  observations  to  failure  were 

t  _  '  *1,20  "  20 '  *2,20  “  X3,20  *  ■Li^’  *4,20  "  0,54  *5,20  “  102  ' 

Estimate  the  aean  life  6  and  give  a  one  and  two-sided  90)1  confidence 

Interval  for  6  based  on  the  data. 

Solution:  this  is  a  non -replacement  situation  with  n  «■  20  and 


r  »  5  •  The  total  observed  life  is  given  by  T. 


\  +  15*5  « 


536  +  2730  *  3266  .  thus  it  follows  fran  (3)  that  '$  •  T^/5  «  3266/5  ■  653  . 
A  two-sided  90)6  confidence  interval  for  6  is  given  by  (357,1657  )  end  a 
one-sided  90)6  confidence  interval  for  0  is  given  by  (409,a> )  .  these 
values  arc  obtained  using  Tables  1  and  2  . 

Exsegils  3:  An  extensive  life  test  has  been  run  a^d  a  0  baaed  on 
r  *<  100  failures  has  been  computed.  Suppose  that  0  •  1000  .  Give  one 


I 
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«nd  two-sided  S%  confidence  intervals  for  9  . 

Solution;  Pros  (10)  the  two-sided  9%  confidence  interval  for  9 

is  given  by 


/  1000 

1000  \ 

1000 

1000  \ 

1  +  h&l 

.  1 ii2£  1  j 

TTI§5 

’  7m 

/too 

Mol 

1 

»  (836,12*4) . 

Proa  (11)  the  one-sided  9%  confidence  interval  for  9  is  given  by 


1  1000  _  t 

/  1000  _\ 

i.idss  ’“I 

*  /r.1555  '  ® 

1  VX55  / 

l  1 

Bxanple  4;  The  total  life  observed  in  obtaining  5  failures  in  9205 
hours.  On  the  basis  of  this  information,  can  we  be  9%  confident  that  the 
probability  of  surviving  for  a  tine  ta  «  100  is  >  .90  1 

Solution:  From  (20)  it  is  known  that  in  order  to  be  99JJ  confident 
that  the  probability  of  surviving  for  a  tins  t*  m  100  is  >.9  ,  it  is 
necessary  that  the  total  observed  life 

T5  -  ^05(10)100/2  log  ^  »  8689  . 

Since  the  total  life  observed  in  obtaining  5  failures  is  9205  hours,  we 
can  answer  in  the  affirmative . 
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i 


FjmupI*  Suppose  that  we  vest  to  keep  a  Mechanism  containing 
lOOb  tubes  in  continuous  operation  for  1000  hours.  Suppose  that  all 
we  know  about  the  tube  lifaiBbaaed-  on  the  data  contained  in  Example  1. 
Based  on  these  data,  how  many  tubes  should  we  expect  to  put  in  as 
replacements  for  those  which  fail  during  the  1000  hour  period?  Find 
a  two* sided  sod  one-sided  95^  confidence  interval  for  the  expected 
nutter  of  replacements  needed. 

Solution:  We  are  in  effect  observing  a  Poisson  process  with  failure 
rate  ■  1000 /B  .  The  maximus  likelihood  estimate  of  \  is,  from 
the  solution  to  (1),  given  by  X  •  1000/0  »  1000/1628  -  6lt  .  Therefor* 
the  expected  nuaber  of  replacements  over  1000  hours  is  given  by  loodl  *>  6lk. 

In  example  (1),  we  computed  (795, 501*0  as  the  two-sided  9 %  confidence 
interval  for  6  .  this  gives  the  two-sided  95^  confidence  interval  for 
the  expected  nunher  of  replacesents: 


I  106  106  ' 

^36K  ’  795 


(199,1258)  - 


In  example  (1),  we  computed  (889,00)  as  the  one-sided  9 5$  confidence  inter¬ 
val  for  0.  Therefore  a  one-sided  95£  confidence  interval  for  the  expected 
nuiber  of  replacements  is  given  by 


0,  $$  J  -  (0,1125)  • 

The  limits  are  very  vide,  because  the  data  are  of  course  very  Inadequate, 
but  they  do  give  us  some  idea  of  what  we  may  expect  to  get. 
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ft— art:  Mon  generally  suppose  wo  want  to  keep  a  eechani—  containing  N 
tubes  la  continuous  operation .  To  do  this  H  good  tubes  aunt  be  in  operation 
at  ail  times.  Suppose  that  we  want  this  condition  to  hold  for  a  tl—  interval 
of  length  T  .  How  eacy  tubes  can  we  expect  to  insert  as  replacements  in  a 
ties  T,  basing  our  estimates  on  one  or  sue  previous  life  testa? 

As  indicated  in  txaaple  5  we  are  in  effect  observing  a  Poisson  process 
with  parameter  K&  *  I/O  .  Therefore,  the  expected  nuaber  of  replacements 
if  we  wish  to  keep  H  ltens  functioning  at  all  tla»3  in  an  interval  of  length  T 
is  given  by  'igT  •  MT/0  .  If  (0^  <  0  <  ®g)  is  a  100(1 -a)  percent  two-sided 
confidence  interval  for  6  ,  then  a  lOO(l-a)  percent  confidence  Interval  for 
the  expected  nuaber  of  replacements  is  given  by  (MT/e^NT/e^)  .  If  (©^,oo ), 

Is  a  lOO(l-a)  percent  one-sided  confidence  interval  for  0  ,  then  a  lOO(l-o) 
percent  confidence  interval  for  the  expected  nuaber  of  replacements  is  given 
by  jo,  • 

'  in  exaaple  5,  a  «  .05,  N  *  1000,  T  •  1000,  #  «  1628,  ei  -  795,62  -  5014, 
and  6^  »  889  . 

Ex— lie  6:  Given  the  data  in  problem  1,  find  a  nuaber  t  such  that  we 
can  assert  with  95$  confidence  that  at  least  90)6  of  the  population  survives  v  . 
(Mote  that  this  is  a  tolerance  statement). 

Solution:  He  noted  in  Remark  (1)  following  our  discussion  of  interval 
estimates  for  the  quantile  xp  that  one-sided  lOO(l-a)  percent  confidence 
statements  regarding  axe  also  tolerance  stata—ots  in  which  we  can  have 
lOO(l-a)  peroent  confidence.  Bence  using  the  solution  to  1(b)  we  can  assert 
that  t  m  93.7  .  Baaed  on  the  data  we  can  assart  with  95f  confidence  that 
at  least  90$  of  the  population  survives  v  •  93.7  hours. 
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To  solve  (b)  graphically  in  the  two-sided  case  we  see  where 
the  horizontal  line  x  »  1000  Beets  the  two  lines  x  «  795  log 
and  9014  log  ^  .  The  two  values  of  p  obtained  are  .28  and 
.82  ,  Hius  a  95$  two-sided  confidence  interval  for  p.  *  . 

»  e'^000/0  (i.e.,  the  probability  of  surviving  1000  hours)  is  given 
by  (.28,  .82)  .  In  the  one-sided  case  the  horizontal  line  x  »  1000 
intersects  the  line  x  «  889  log  i  at  p  *  .32  .  Hence  we  can 
state  that  (.32,1)  is  6  95$  confidence  interval  for  p^»  iqqq  ■ 

»  e-1000/®  .  '  *  1  ^ 


HlO. 


'i 


Section  2. 

Aa  Estimation  Problem  (Fixed  time  of  truncation.  I  tew  which  fail  are 
replaced  by  new  item-) 

Problem:  n  items  are  placed  on  life  test  at  time  t  «  0  .  As  the  test 
proceeds,  items  which  fail  are  replaced  by  new  items.  Life  testing  is  terminated 
at  time  t*  .  It  is  assumed  that  the  underlying  p.d.f.  of  life  is  given  by 

f(tjfl)  -  |  ,  t  >  o  ,  e  >  o  . 


We  wish  to  do  the  following: 

(I)  Estimate  0  . 

(ii)  Make  one  and  two-sided  confidence  statements  about  0  . 

(ill)  Make  probability  statements  about  the  proportion  of  items  having 
life  greater  than  t*  . 

Solution:  In  what  follows  let  r  *  number  of  items  which  fail  in  (0,t#), 
then  the  solutions  are  as  follows: 

(i)  The  maximum  likelihood  estimate  for  6  is  given  by  at*/?  * 

(II)  A  one-sided  lOO(l-a)  percent  confidence  interval  for  0  is  given 


by 

(1) 


gnt 

^(2r  +  2)  ’ 


ao 


\  a 


A  two-sided  100(1 -a)  percent  confidence  interval  for  9  is  given  by 


(ill)  From  the  results  in  (il)  regarding  the  one-sided  lOO(l-a) 
percent  confidence  intervals  for  6  we  can  be  lOO(l-a)  percent  confident  that  at- 


-3- 18- 


least  luO  bjl  of  the  population  survives  t*  hours,  with 


bee 


(2r+2)/2n 


In  other  words  a  lOO(l-a)  percent  one-sided  confidence  interval  for 


b*e 


is  given  by 


(2r+2)/2n  \ 

e  >  1  I  • 


Fran  the  results  in  (ii)  regarding  the  two-sided  iOO(l-a)  percent 
confidence  intervals  for  8  ,  ve  can  say  that  if  ve  observe  r  failures 
in  (0,t  )  then  a  two-sided  lOO(l-a)  percent  confidence  interval  for 

)j  K  tt  'll*  nlmn  hv 


is  given  by 


->?  (2r+2)/2n  -X8  0  (2r)/2n 

•  5  .  .  l'5 


Proof*  Essentially  ve  are  observing  a  Poisson  process  with  parameter 
X*  a  nx,  -  where  X  »  ^  -  If  we  observe  r  failures  in  (0,t  )  then  the 
moviana  likelihood  estlsnte  for  is  given  by 


r 

■  t» 


V  r 
*  *  5  “5T 


(7) 


Therefore 


(8) 


*.±. 


at 

r 


and  this  establishes  (i)  . 

It  can  be  shown  that  ths  probability  of  observing  r  or  fewer 
failures  in  (0,t*)  ,«e  given  by 


(9) 


Pr(k  <  r|0)  »  _> _ e  ^(nt*/©)Vk  2 

fc*0 


Thus,  if  0  <  2nt*/X^(2r+2)  then  Pr(k  <  r|fl)  <  a  .  This  inplies 
that  if  v*  observe  the  event  k  a  r,  then  we  are  lOO(l-a)  percent 
confident  of  the  correctness  of  the  assertion  that  8  >  2nt  /  a ,  (2r+2). 

In  a  slnilar  way  it  can  be  shown  that  if  8  <  2nt*/X?(2r+2)  than 

2 

Pr(k<  r|0)<  |  and  if  Q  >  2nt*/*2  a(2r)  then  Pr(k  >  r|0)  <  |. 

1-  2 

From  this  it  follows  that  if  we  observe  the  event  k  «  r,  then  ve  are 
lOO(l-cr)  percent  confident  of  the  correctness  of  the  assertion  that 


Remark  1:  Define  9  ac  6*8  «  nt*/i*-!-l  .  Then  one  can 

write  the  one-sided  lOO(l-a)  percent  confidence  interval  for  6  as 

9  t  as  I  and  the  two-sided  lOO(l-a)  percent  confidence  intenml 
,<[2r+2] 


A 

2  r  © 


12**1  8  _ _ _ 

^(2r+2)  r  j2r)l 


!.  2' 
A  2 


Urns  9  is  involved  in  computing  the  one-sic_d  Interval  and  in 

A 

finding  the  left-hand  end  point  of  the  two-sided  interval.  9  1b 
Involved  in  finding  the  right-hand  end-point  of  the  two-sided 
interval.  It  la  now  clear  that  we  can  use  Tables  1  and  2  in  order 
to  compute  the  confidence  intervals. 

Remark  2:  If  r  *  0,  only  the  estimator  0  makes  sense  and 
only  one-sided  intervals  of  the  fona  (1)  should  be  used. 

Remark  3:  The  two-sided  confidence  intervale  ifcr  9  given  toy 
formula  (2)  are  direct  cons  a  queues  of  formulae  for  two-sided  con¬ 
fidence  intervals  for  the  parameter  \  in  a  Poisson  process  given 
by  P.  Garwood  in  Biometrlka  28,  437-4*2,  1936,  This  question  is 
also  treated  in  S.8.  Pearson  and  H.O.  Hartley,  Bloastrika  Tables 
for  Statisticians,  Vol.l,  pp. 74-77,  Cartridge  itoiversity  Praaa,1954. 


(13) 


( 
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2nt*  Ior  -  2nt*  Ion  i 

)g(2r+2)  X2  a(2r) 

5?  I" 


-  ^  (ggglllog  |  .  2r£log  | 


>jS+2) 

5 


p  # 


v(2r) 


Fonsulse  (12)  and  (13)  given  lOO(l-a)  percent  one  and  two* sided 
confidence  beads  for  the  entire  distribution. 

Beaark  7:  It  follows  froa  (12)  that  if 


w 


_  2nt  1 

t  a  - -  lo*  - 

*J(2r+2)  p 


then  we  can  assert  with  lOO(l-c)  percent  confidence  that  (i,oo )  is 
a  100  p  percent  tolerance  interval.  More  precisely,  if  one  observes 
r  failures  in  (o,t*)  (where  n  item  are  constantly  kept  on  test) 
then  we  can  be  lOO(l-a)  percent  confident  that  the  probability  of 
surviving  for  at  least  tine  t  is  >  p  (or  that  the  fraction  of 
the  population  surviving  t  or  aore  hours  is  >  p).  in  terns  of 
e,  (12)  can  also  be  written  as 


(15) 


T  x 


(fj+SiJ!  log  i  . 

^(2r+2)  P 


this  again 


it  eeay  to  ccapute  t  tram  Tables  2  and  3. 


1  — 
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It— ark  8:  lb  the  section  devoted  to  the  testing  of  hypotheses 

we  stalled  a  truncated  replacement  life  test  procedure  of  the  following 

kind,  n  item  are  placed  on  test/  end  it  Is  decided  la  advance  that 

the  experiment  will  he  terminated  at  min  (t  :t*),  where  v  Is 

r0/B  ro,n 

a  rand—  variable  equal  to  the  tins  at  which  the  rj *h  failure  occurs 
and  t*  is  a  truncation  tine,  beyond  which  the  experiment  will  not  be 
run.  Both  rQ  and  t*  era  assigned  la  advance  before  life  testing 
starts.  If  the  experiment  Is  terminated  at  t  „  (i.e. ,  if  r 

TmU  C 

o 

failures  occur  before  time  t*),  then  the  action  in  terms  of  hypothesis 
testing  is  the  rejection  of  so—  specified  null -hypothesis .  if,  however, 
the  experiment  is  terminated  at  time  t*  (i.e.,  the  r  **  failure  does 
not  occur  before  time  t*),  then  the  action  In  teres  of  hypothesis 

•an* 

testing  is  the  acceptance  of  some  specified  null-hypothesis. 

8vgpose  now  that  such  a  test  has  been  run  and  that  we  would  like  to 
see  the  data  obtained  not  only  for  testing,  but  also  for  estimation.  It 
la  g—arslljr  recognised  that  there  are  difficulties  associated  with 
using  such  data,  since  the  stopping  rule  usually  affects  the  estimates 
“hi oh  o—  be  obtained.  It  is  interesting  to  point  out  that  for  the 
truncated  life  teat  under  discussion  the  following  rule  gives  lOO(l-a) 
percent  one-sided  confidence  intervals: 

U)  I*  T_  >  t#,  i.e..  If  the  mnfcer  of  observed  failures  k  In 
«  ° 

(0,t  )  is  0*1,2,  ....  ,re  -1,  then  a  one-sided  100(1 -a)  percent  con- 
fid— os  Interval  is  given  by 


i 


(16) 
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20^ 


xJ(2k+2) 


(ii)  If  Tr  <  t*  ,  a  one-sided  lOO(l-a)  percent  confidence  interval 
o 

is  given  by 


(17) 


;  2nt 


r 

o 


*f(2r  ) 

a'  o 


In  Appendix  3  E  we  prove  that  equations  (16)  and  (17)  generate 
100(1-0 )  percent  one -aided  confidence  Intervals. 

One  night  conjecture  that  two-sided  100(1 -o')  percent  confidence 
Intervals  can  be  defined  In  an  analogous  way  as: 


(16) 


and 

(19) 


i  2nt  \ 

|  "y*-""  ,  a»  \  if  k  c  0 

(  ?<2>  ) 


k  -  .1,2, 


2nx 


X  a(2r  > 

1-  |  ro  l 


if 


Tr  5  t*  . 


We  have,  up  to  now,  not  been  able  to  establish  this  conjecture  rigorously. 
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numerical  Examples 

Problag,  1;  30  Item  are  planed  on  teat.  Item  which  fail  are 

replaced .  'Ehe  life  test  is  shopped  after  100  hours  have  elapsed. 

Five  failures  are  observed  in  the  course  of  the  experiment .  Assigning 
that  the  underlying  distribution  of  life  is  exponential,  find 

(a)  An  estimate  of  the  mean  life  8  .  Giva  one  and  two-sided  99$ 
confidence  internals  for  6  . 


(b)  Make  one  and  two-sided  95)6  confidence  statements  for  the 
probability  of  surviving  100  horns. 

(c)  Make  one  and  two-sided  95$  confidence  stt.tejueats  for  the 
probability  of  surviving  50  hours, 

Solution: 

(a)  In  this  problem  n  »  30,  t*  *=  100,  the  observed  number  of 

failures  is  r  «  5.  Thus  the  maximum  likelihood  estimate  for  8  is 
given  by  8  *•  nt‘  /r  «  3CCO/5  <=  6CC.  GubstlvUtxQg  in  formula  (i)  and 
using  (12)  *  21.026,  one  gets  the  one-sided. 95$  confidence 

interval  (285,  00)  .  Substituting  in  formula  (2)  and  using  >?Q 2j(12)«23.337 
.24?  one  gets  a  95$  two-sided  confidence  interval  (257,1048). 

(b)  A  one-nided  95$  confidence  interval  for  surviving  t#  «  100 
hours  is  given  by  (e*21*026/60  ,1)  «  (e”3504  ,1)  «  (. 704,1). 

A  two-sided  95$  confidence  interval  for  surviving  t.*  «  xoo  hours 
ir,  given  by 

(„-23.33T/60_  ,-3.247/60,  .  (.-3889,  ,-0541,  ,  ^  _ 

(c)  One  end  two-sided  95$  confidence  intervals  for  the  probability 
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of  surviving  t  -  50  hours  are  given  by  (e  »  (. 8393*1)  ao& 

(e-.19^5}  c-*0271)  m  (.8232,  .9733)  respectively. 

Problea  2:  Given  the  data  in  example  1.  Ectimate  t  so  that  we 
will  be  95$  confident  that  the  probability  of  surviving  t  hours  is 
at  least  .9  . 

Solution: 

a  •  30,  t*  »  100,  r  «  5,  a  «  .05,  and  p  *  .9  .  Thus  substituting 
in  (lb)  we  have 

’  "  I®  -  30-1- 

More  directly,  using  tables  2  and  3  and  noting  that  9  »  6  =  -500, 

one  gets  1  •  (500)( -571)( -lC^)  *  30-1  . 

On  the  basis  of  the  data  we  can  be  95$  confident  that  the  probability  of 
surviving  t  -  30.I  hours  is  >  .9  . 

Problea  3:  A  truncated  replacement  test  consists  of  placing  30  items 
on  test  for  at  aost  100  hours.  If  3  failures  occur  before  100  hours,  the 
life  test  is  stopped  at  once  and  the  lot  is  rejected.  If,  however,  3  Items 
have  not  yet  failed  by  the  time  100  hours  have  elapsed,  the  teat  is  termi¬ 
nated  at  100  hours  with  acceptance.  Items  which  fail  axe  replaced  at  once 
by  new  iteas.  Give  a  95$  one-sided  confidence  interval  for  G  if  one 
observes  exactly  one  failure. 

Solution:  We  use  formula  (16)  with  It  »  1,  hence  a  one-sided  95$ 
confidence  interval  is  given  by 


t 

1 


r 
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m  {632  j  GO  )  - 


Problot  kt  Suppose  It  happened  that  the  third  failure  was 
observed  to  occur  at  50  hours .  (Jive  9 9#  one  end  two-sided  confidence 
Intervals  in  this  case. 

Solution:  WO  use  foxmla  (17)  with  k  «  3t  *  50  .  >' 

"  ^  ,3  9  jw 

Hence  a  one-sided  95#  confidence  Interval  is  given  by 


( 2nt_ 


OO  I  »•» 


00  |  «  (238*  00 } 


Substituting  In  fonula  (19)  ve  get 


( .  fvl  - !  sv  ’  asr-V  <■**«» 

\  .025<6)  **975  R  ,M  / 


as  a  two-oidcd  95#  confidence  interval. 
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SECTION  3 

AS  ESTIMATION  PROBLEM 

(Fixed  time  of  truncation  t*;  failed  itene  are  not  replaced) 

Problem:  n  items  are  placed  on  life  te>*t  for  a  time  t*.  At 

—  . ■“ 

the  end  of  this  time  one  counts  the  number  of  items  that  have  failed  in 
the  time  Interval  £o,  t*J.  Items  that  fail  are  not  replaced.  We  wish 
to  do  the  following; 

(i)  Give  an  estimate  for  the  probability  of  surviving  for  a 
length  of  time  t*  and  further  estimate  the  mean  life  9,  if  the 
underlying  distribution  is  exponential. 

(ii)  Make  one  and  two-sided  confidence  statements  about  the 
probability  of  living  for  more  than  t* .  Stated  in  reliability  language 
we  wish  to  make  probability  statements  about  the  reliability  of  items 

in  [o,  t*J, 

(ill)  Make  one  and  two-sided  confidence  Bt&temects  about  the  mean 
life  9  in  the  case  where  t.he  underlying  distribution  is  exponential. 

Solution:  In  what  follows  let  r  =  number  of  items  which  fail  in 
to,  t*J,  then  the  solutions  are  as  follows: 

(i)  The  maximum  likelihood  estimate  of  the  probability  of  sur¬ 
viving  more  than  t*  time  units  is  given  by 

*  A  -t*/@ 

If  the  underlying  distribution  is  exponential,  then  p  »  e  '  and  hence 


(2) 


.  ■e*e*nL&mt*mr**j*m*  v. 


X 
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(ii)  There  is  a  confidence  of  100(1  -  c£  )  percent  attached  to 
the  statement  that  at  least  100  b$  of  the  population  survives  for  a 
length  of  time  t*  with  b  given  by  the  formula 


(3) 


_ 1_ _ 

1  +  (Hr)  p<* (2r  +  2> 20  * 2r) 


In  other  words  the  one-sided  100(l  -  eX )  percent  confidence  Interval  for 
the  probability  of  surviving  t*  time  units  Is  given  by 


(nl>  “2^  18  deflned  in  a  that  ^(^(n^,  ng)  2  F^Cn^,  ng))  »  ct, 

where  ng)  is  the  F  distribution  with  n^  degrees  of  freedom  in 

the  numerator  and  ng  degrees  of  freedom  in  the  denominator. 

A  two-sided  100(l  -  at)  percent  confidence  interval  for  the 
probability  of  surviving  t*  time  units  is  given  by 


(4) 


(l*(r r?)V^2- 


a>-2r)  Wrrfn)  p,o.(2r-  a-2"s) 


>4 


;> 


These  results  sre  completely  distribution  free. 

(ill)  In  the  case  where  the  underlying  distribution  is  exponential. 


one  obtains 


as  a  one-sided  100(1  -  <£)  percent  confidence  interval  for  0  and 


n 


<!► 
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(6) 


t* 


F  (2r+2,  2n-2r )J  log  £  1  +^f~TT>)  Fi^2r^2n"2r+2J 


as  a  two-sided  100(1  -  oL)  percent  confidence  interval  for  0. 

Proof i  (i)  Formulae  (l)  and  (2)  are  obvious. 

(ii)  Suppose  that  we  observe  r  failures  in  the  time 
interval  [o,  t*J  and  that  p  *  probability  of  failing  in  £o,  t*J  and 
q  -  1  -  p  =  probability  of  surviving  in  [o,  t*J.  Suppose  that  pQ  is 
such  that 


(7) 


r~»  / u>  k  n-k  , 

L  (k>  pD  %  “ 

k*o 


then  we  can  state  that  Pr(k  <  r  /  p)  £  at.  if  pipQ.  Hence  if  we  observe 
k  =>  r,  we  can  be  100(1  -  d)  percent  confident  that  p  <  pQ  or  that 
q  a  1  -  p  >  1  -  Pc  =  cjQ.  The  question  arises  as  to  how  one  computes  qQ. 
This  can  be  done  very  easily  by  expressing  (7)  ati  an  incomplete  Beta  Func¬ 
tion  and  then  using  the  well-known  relationship  between  the  Beta  and  F 
distributions.  If  this  is  done,  one  discovers  that 


i 


(8) 


% 


1  +  ir^r)  +  2,  2n  -  2r)  ' 


where  F^  (n^  n2)  is  defined  in  such  a  way  that  Pr^F^,  ng)  >  (n^, 

*  ot,  and  where  n^  and  cure  the  number  of  degrees  of  freedom  in  the 
numerator  and  denominator  respectively.  Thus  (3)  is  established.  In  this 
connection  one  should  also  read  S.  Takada  and  S.  Shimada,  Bart  1,  July  195^, 
pp„  147  and  151*  See  bibliography  given  in  the  Appendix, 

In  Table  4  we  give  the  values  of  qo  for 
n  =  1(1)20(5)30(10)50(25)100(50)200(100)500;  for  oL  =  .01,  .05,  .10,  .25,  .50 


I 


4 


I 


i  Kr 


tv- 

«► 
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and  r  *  l(l)  .  min  (n.  20) „  In  Table  5  we  tabulate  q„  for  n  =  1000, 

o 

5000,  10000,  50000,  100000,  500000,  1000000;  for  oL  =  ,01,  ,0?,  .10, 
.25,  .50  and  r  =  l(l)20(l0)l00,  200,  500. 

Two-sided  confidence  results  ore  obtained  by  finding  p  =  and 
p  =  Pp  such  that, 

k  n-k 


(9) 


and 


(10) 


£f?XX\  it  t 

(k)  Pi  qx 

k*o 


oL 

2 


n 


E  ^  p2  q2~  -  2 

k*r 


Hence  if  k  *  r  is  observed  we  can  be  100(1  -  cL)  percent  confident  that 
P2  <  P  <  Px  or  that  qx  <  q  <  q£. 

The  computation  of  q1  and  q2  involves  expressing  (9)  and  (10)  as  an 
incomplete  Beta  Function  and  then  using  the  well-known  relationship  be¬ 
tween  the  Beta  and  F  distributions.  If  this  is  done,  it  turns  out  that 


(11) 


i  +  <■ 


2n  -  2r) 


and 


1  + 


( n  -  r  +  l)  T1j*  *2r>  2n  -  2r  +  2) 


Thus  (4)  is  established. 

Tables  for  and  are  being  computed  for  the  values  of  n, 
el,  r  used  in  Tables  4  and  5* 


* 

r 


“332- 


lu  the  particular  case  where  the  underlying  distribution  happens 
to  be  exponential,  (6)  implies  that  (5)  is  a  one-sided  100(l  -  06)  percent 
confidence  interval  for  9  and  (ll)  implies  that  (6)  is  a  two-sided 
100(l  -  ot)  percent  confidence  interval  for  0. 

Remark:  One  end  two-sided  100(1  -  6L)  percent  confidence  intervals 
for  the  probability  of  surviving  an  arbitrary  time  T  not  necessarily  * 
t*  are  given,  in  the  exponential  case,  by 


(12) 


(2r  +  2, 


2n  -  2r)J 


and 


<*)  N1  +  (Hr)  y#(2r+2>  2n_2rj  ^ '  [x  fi^>  *-*>*>J 


respectively. 

Remark;  It  happens  sometimes  that  n  is  very  large  and  r  is 
very  small.  It  is  useful  t~>  note  that-  in  this  case  (\)  becomes 


b 


1+  (S-LJ.J  F*(2r  +  2,  co) 


In  other  words,  the  one-sided  100(1  -  cL)  percent  confidence  interval  for 
the  probability  of  surviving  time  t*  is  given  by 


1  _ 

1  +  "Xa(?£.+  2) 

2d 


Similarly  a  two-sided  100(1  -  ot)  percent  confidence  interval  for  the 
probability  of  survival  is  given  by 
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2 

In  (15)  and  (16)  we  use  the  fact  that  (ra,  oo)  *  X 

In  all  of  the  results  obtained  up  to  this  point  in  this  section, 
we  have  not  made  any  use  of  the  failure  tines  of  those  items  which  did 
indeed  fail.  Because  of  this  we  were  able  to  state  a  certain  number  of 
non -parametric  results.  However,  in  the  event  that  the  underlying  dis¬ 
tribution  of  life  really  io  exponential  we  are  clearly  losing  some  informa¬ 
tion,  at  least  when  n,  the  number  of  items  originally  placed  on  test  is 
small.  We  say  this  because  of  the  fact  that  If  n  were  very  large,  then 
we  would  be  effectively  dealing  with  a  replacement  situation.  In  this 
case,  the  knowledge  of  actual  times  to  failure  is  Irrelevant  if  the 
underlying  distribution  is  really  exponential.  It  is  only  for  small  or 
moderate  sizes  of  n  that  it  would  make  a  difference  whether  or  not  we 
use  our  knowledge  of  the  actual  times  to  failure  of  those  items  which 
did  fail. 

Throughout  we  assume,  as  before,  that  we  start  the  life  test  with 
n  items  and  that  we  do  not  replace  failed  items.  Let  r  =  number  of 
items  which  fail  in  (0,  t*)  and  let  ^  be  the 

failure  times.  We  assume  further  that  the  underlying  distribution  is 
exponential.  An  exact  solution  to  the  problem  of  finding  100(1  -  M.) 
percent  confidence  intervals  for  6  is  easy  in  principle,  but  difficult 
to  carry  out.  Hence  we  give,  without  proof,  some  approximate  procedures 


»  -  'JT' 


*  • 


i 

I 

I 


<v 

O 
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which,  are  good  enough  in  most  practical  problems.  The  first  result  is 
that  approximate  one-sided  100(1  -  <J.)  percent  confidence  intervals  for 
9  are  given  by 


(17) 


/.  3(f)  ,m) 
V  ^(2r  .  2)  / 


for  r  s  0,  1,  2,  . n  ~  1 


A 

where  T(t*)  *  2j  +  (n  “  r)t* 
i«l 


and 


*l(a>) 


,  00 


where 


Krn)  -  £  7;  ,  if  r  «  n, 
i«l 

Approximate  two-sided  100 (l  -  eL)  percent  confidence  intervals  for 
9  are  given  by 

/ 


(18) 

by 


(  -2S2L  ,  CD  ) 


if  r  *  0 


2T(t») 


2T(t*) 


*X|(2r  +  2)  7q_|(2r) 


;) 


if  r  *  1,  2,  n  -  1 


and  by 


a(rn; 


if  r  «  n. 


Formulae  (17)  and  (18)  should  be  compared  with  (5)  and  (6)  respectively. 

Remark:  It  is  convenient  to  define  6*  as  Q  *  9—  * -»  Kill , 

r  +  1  r  +  1  1 

for  r  *  0,  1,  2,  . n  -  1  and  as  0  »  0  =  T(?'n)/n  for  r  =  n.  Formula 
(l?)  then  becomes 


!»w<ggggiigjgjl BagjjPS 


(17') 
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,  oo  l  for  r  «  0,  1,  2, 

X^(2r  +  2)  J 


n  »  1  and 


i  ^y  foi 


sr  r  =  n 


and  formula  (18)  becomes 


(18') 


— . .  ,  *  )  for  r  «  0 

*>  J 


*  2}  >Cf_*(ar) 

V.  2 


for  r  -  1,  2,  . . . ,  n  -  1 


’X^&O  Xl-< *(2n)y 

2  2  / 


for  r  -  n. 


As  was  done  before,  let  us  define  the.-  quantile  an  that  life 
such  that  Pr(X  a  x  )  »  p,  i.e.,  x  s  a  log  ~  Then  approximate  one- 
and  two-sided  100(1  -  ot)  percent  confidence  intervale  for  x  are  obtained 
by  multiplying  the  formulae  in  (17),  (17*),  (18),  (l8‘)  by  log  ~  .  Further¬ 
more  it  follows  that  if 

/V 

(19)  t  -  -l|OiU-e-.  log  i  ,  r  =  0,  1,  2,  n  -  1 

P  Aoc/2r  +2)  p 

and  A 

=  -7^  -  log  i  for  r  -  n 

Xl<2n)  * 

then  wo  can  be  approximately  100 (l  -  ot)  percent  confident  of  tie  correct¬ 


ness  of  the  assertion  that  the  fraction  of  the  population  surviving  t  or 

Jr 
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t 


more  hours  is  p.  As  befoi-e,  one-sided  confidence  intervals  on  quantiles 
are  eevi. . 1  lent  to  one-sided  tolerance  statements  about  tbe  population, 
with  the  same  confidence. 

•7  1 

If  we  want  to  place  approximate  one  and  two-sided  100(1  -  «L) 

■  •*  T/e 

percent  confidence  intervals  on  •  e  '  ,  the  probability  of  sur¬ 

viving  TT  hours,  the  results  are  as  follows : 

TXl(2rr2)  \ 


(20) 


and 


* *  *  , 1)  0,  f  ) 


2)/2T(t*) 

) 

for  v  *  0,  1,  2,  n  -  1 


-  ?*;j2n)/2T(Tn) 


'  7  1 


or  r  =  n 


are  approximate  100(1  -  at)  percent  one-sided  confidence  intervals  on 

T 


and 


(21) 


life) 

2nt*  1 
e  ,1 


for  r  »  0 


and 


1 

\ 


l 


7-X^(2r  +  2) 


2{r  +1)?" 

e  > 


t 


rXl-et  (2r) 

2 


T7CV  <jt(2n) 
2 


for  r  «  1,2, 


for  r  a  n 


2il0-  2rf9 

»  e 

are  approximate  1C0(1  -  cL)  percent  two-sided  conf  idence  internals  on  px  „ 


m&nm 


OO*  *■ 


Suppose  that  the  data  are  originally  obtained  in  the  course  of 

,  ! 

running  a  truncated  non-replacement,  life  teat  wit.n  preassigucdi  truncation 
time  t*  and  maximum  allowable  number  of  failures  r  .  The  stopping 


is  a  random  variable  equal  to  the 


rule  is  min(T  t*)  where  T 

ro-n'  ,  ro,n 

time  at  which  the  rQ  *t.h  failure  occurs .  Then  the  following  rule  gives 

approximate  100(l  -  cL)  percent  confidence  intervals  for  8. 


(i)  If 


T  >  t*,  i,e,,  if  the  n!jmber  of  observed  failures  k 


in  (0,  t*)  is  0,  1,  2.  . r  -  1  then  an  approximate  one-sided 
100(1  -  oL )  percent  confidence  interval  is  given  by 


(22) 


wnere 


(  2T(t») 


— __  >  co 


^  X^(2k  +  2) 


) 


0.  1,  2., 


x. 

T(t*)  =  y  r.  +  {n  -  k)t* 

'  fci  A 


r  -  1 
o 


or  equivalently  as 


t(k  v 

Ki2k  * 2) 

where  «  «  T(t*)/k  +  1. 

(ii)  If  T  t*,  then  the  appropriate  interval  is 


(23) 


2T(tt  ) 


A 

2r  C 
o 

xa<2ro> 


00 


where 


T(rr  ) 

o 


- 1 


i»i 


Ti  +  ’n  ■  toH* 


and  6  c  TfTr  )/rQ 
*0 
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In  an  analogous  way,  we  can  obtain  approximate  two -Bided  100(1  -  a) 

Percent  confidence  intervale  for  6. 

Remark:  We  wish  to  reemphasize  that  in  the  last  few  pages  ve 
hare  given  results  which  have  not  been  and  probably  cannot  be  rigorously 
established.  However,  they  can  be  used  as  good  approximations  to  true 
results.  Further  discussion  of  this  point  is  deferred  to  the  Appendix. 

Numerical  Examples 

Problem  1:  20  items  are  placed  on  life  test  for  100  hours.  Two 

items  fail  before  this  time.  Items  which  fail  are  not  replaced. 

(a)  Make  non-parametric  9 5$  confidence  statements  (one  and  two- 
sided)  about  the  probability  of  items  surviving  100  hours. 

(b)  If  the  underlying  distribution  is  exponential  find  one  and 
two-aided  95$  confidence  Intervals  for  8,  the  mean  life. 

(c)  If  the  underlying  distribution  is  exponential,  give  one  and 
two-sided  95$  confidence  intervals  for  surviving  t  -  50  hours. 

Solution:  (a)  In  thi3  problem  n  *  20,  r  ■*  2,  a  s*.  .05,  t*  *  100.. 
Since  t  q^(6,  36)  •  2  > 36 j  it  follows  from  (33 )  that  a  one-sided  95$ 
confidence  interval  for  the  probability  of  surviving  t*  »  100  hours 
is  given  by  (.718,  l).  Since  F ,025(6*  36)  »  2.79  ®nd 
F  2^(4,  36)  *  1/9.42,  it  follows  from  (4)  that  a  two-sided  95$ 
confidence  interval  for  the  probability  of  surviving  t.Q  -  100  hours 
la  given  by  (.683,  .986); 

(b)  From  (5)  a  one-sided  95$  confidence  interval  for  9  is  given 
by  (302,  00 )  and  from  (6)  a  two-aided  95$  confidence  interval  is 
given  by  (262,  8C5). 

(c)  From  (12)  a  one-sided  95$  confidence  interval  for  the  probability 
of  surviving  1  »  50  hours  la  given  by  (.6447,  i)  and  from  (13)  a  two-sided 
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95 5  confidence  Interval  for  surviving  T  «■  50  hours  Is  given  by  (.826,  .99*0* 
Problem  2:  Out  of  10,000  items  tested,  no  items  were  observed  to 
fall.  Give  a  one-sided  955  confidence  interval  for  the  probability  of 
survival. 

Solution:  In  this  case  n  *  10,000,  r  =  0.  Since  n  is  very 
large>  F  Q^(2,  20000}  ~  F^(2,  00)  *  3.00>«a4  eo  the  cna«sidod  9*5  con* 
fidence  interval  is  given  by  (.9997,  l).  In  other  words,  we  have  955 
confidence  in  the  assertion  that  the  true  probability  of  survival  1b 
2:  .9997,  if  no  items  are  observed  to  fail  in  a  seunple  of  10,000.  The 
answer  can  also  be  found  very  easily  by  using  Table  5- 

Problem  3:  Out  of  10,000  items  tested,  10  items  were  observed  to 
fail.  Give  one  and  two-sided  955  confidence  intervale  for  the  probability 
of  survival. 

Solution:  In  this  case  n  =  10, COO,  r  *  10.  Since  n  is  very 
large  F  Q5<22,  19960)  ~  F^(22,  ao)  =  1.5^,  and  so  the  one-sided  955 
confidence  interval  for  the  probability  of  survival  is  (.9983,  l) *  In 
other  words,  we  have  955  confidence  in  the  assertion  that  the  true  proba¬ 
bility  of  survival  is  2  .9983,  if  ten  items  are  observed  to  fail  in  a 
sample  of  10,000- 

In  the  two-sided  case  F>025(22,  19§8o)=  I.67  and  F  ^(20,  19982) 
m  1/2.00  »  .500  and  so  the  two-sided  955  confidence  interval  for  the 
probability  of  survival  is  given  by  (.9982,  .9995)- 

Problem  k:  A  sample  of  20  tubes  is  placed  on  test.  Experimentation 
is  truncated  at  time  t*  »  $00 .  Items  which  fall  are  not  replaced.  In 
this  particular  sample  6  items  fail  before  t*  «  500  hours.  The  total 
life  of  the  6  Items  which  failed  before  t*  »  500  was  956  hours.  Estimate 
the  mean  life  3  and  give  one  and  two-sided  955  confidence  statements  for 


6,  the  mean  life 
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Solution:  (l)  Let  us  solve  the  problem  ignoring  the  information 
that  the  total  life  of  the  6  failed  items  *  95^-  Thus  wo  use  formula 
(5)  with  t*  =  5Q0,  u  «  20,  r  =  6,  CL  =  .05.  This  gives  the.  one-sided 
95%  confidence  interval 


CD  = , 


log  ^1  +  ^  F>05(14,  28 )]'  j  (log  ^1  +  |{2.06)J 

&y  -  “)  ■  (■$&  ’  “J  '  <706' 


,  00 


log 


Similarly  using  formula  (6)  we  get  the  two-sided  95%  confidence  interval 


_ _ _ _ .) 

108  (  1  +  A"  F.025^1  28^J  103  (  1  *  Tjf  F.97^12j  30)]  j 


500 
JK6** 


JS22_ 


500 


Vlos  {  1  *  l(2*:i7))  log  f1  2^5}j- 


500  500 

log  (2.1157  *  log'fiTissT, 


m-  m)  •■<*> ■ 


3950). 


(2)  If  we  use  the  fact  that  the  total  life  of  the  6  observed 
failures  ■  956  we  can,  use  (17)  or  (17')  to  find  a  one-sided  95%  confidence 
interval.  In  this  problem 

6 

x(t*>  =  £3  t:  +  i4t*  «=  956  +  7000  .  7956. 

i=rl 

Further  8  *  =  1137  Using  (17')  and  Table  2  we  get  the  one-sided 

95%  confidence  interval  ((1137) (-591),  00)  *  (672.  ®).  Substituting 
in  (18 !)  and  using  Table  1  we  get  the  two-sided  95%  confidence  interval 
^ ( 1137 K  536),  (1326) (2.725)]  *  (609,  3613)-  The  confidence  intervals 


k' 
*  - 

If 
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obtained  by  the  two  methods,  one  of  which  ignores  the  actual  failure 

times  are  surprisingly  close  considering  the  smallness  of  the  sample 

and  the  fact  that  we  are  dealing  with  a  specific  experiment. 

Problem  5’  Given  the  data  of  Problem  5,  find  one  and  two-sided 

95$  confidence  intervals  for  x  Q,  the  time  which  is  such  that  90  per- 

*y 

cent  of  the  items  in  the  population  live  longer  than  x 

Solution:  We  multiply  the  numbers  obtained  as  95$  confidence 
limits  for  9  by  log  Thun  the  one-sided  9556  confidence  interval 
for  x  Q  is  given  by  (71,  oo  }  and  the  two-sided  95$  confidence  inter- 
val  is  given  by  (64,  381) * 

Remark:  We  can  interpret  the  one-sided  9556  confidence  interval 
for  x  ^  as  a  one-sided  tolerance  statement,  namely  on  the  basis  of  the 
data  we  can  be  95$  confident  in  making  the  assertion  that  at  least  90^ 
of  the  population  survives  x  Q  ■  71  hours. 

Problem  6:  A  certain  company  guarantees  a  television  tube  for  the 
first  month  of  use.  Out  of  1000  tubes  sold,  50  are  returned  under  this 
guarantee . 

(l)  Make  a  non-parametrlc  one  and  two-sided  confidence  statement 
about  the  proportion  of  tubes  lasting  at  least  one  month. 

(ii)  Assuming  the  exponential  distribution  to  be  valid,  estimate 
the  mean  life  0,  and  give  one  and  two-sided  95$  confidence  intervals 
for  0 .  | 

(ill)  Assuming  the  exponential  distribution  to  be  valid,  estimate 
x  g  ,  the  time  when  we  may  expect  50$  of  the  tubes  to  have  failed.  Flace 
one  and  two-sided  95$  confidence  intervals  on  x  c' 

Solution:  Clearly  this  problem  can  be  considered  as  a  truncated 
without  replacement  situation  with  n  *  1000,  t*  =  1,  and  r  »  50.  The 


L 
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problem  can  also  be  considered  as  consisting  of  1000  non-replacement 
truncated  life  tests,  where  each  life  test  consists  of  testing  n  =  1 
item  for  at  most  t*  «  1  hour.  The  customer  carries  out  this  life  test 
and  in  a  sense  accepts  (keeps)  the  tube  if  it  survives  for  one  month 
and  rejects  the  tube  (is  given  a  new  tube)  if  the  failure  occurs  before 
one  month.  We  will  assume  that  accurate  records  have  not  been  kept  and 
that  we  must  base  our  estimate  on  the  number  of  failures  reported .  From 
(i)  the  maximum  likelihood  estimate  of  the  proportion  of  tubes  surviving 
t*  »  1  month  is  given  by  'p  »  ~  *950.  A  one-sided 

95^  confidence  interval  for  the  probability  of  surviving  t*  ~  1  month 
is  given  by  substituting  in  (3  * )  with  r  -  50,  n  =  1000.  This  confidence 


interval  is 


1  *  95§  F.05<102' 


1  - 


1  +  (1.25) 


,  1  J-  (-937,  1). 


A  two-sided  95$  confidence  interval  is  given  by  substituting  in  (9).  This 


gives 


+  950  f.025(102'  1900)  1  +  oir  F  W100»  1902) 


951  .975 


1  +  — Si  (l«3l)  ’  1  +  —52  — i. 

950  951  1.35 


(•93^,  -963). 


This  gives  the  solution  to  (i). 

To  solve  (ii)  we  substitute  in  (5)  and  (6)  respectively.  This  gives 
the  one-sided  95^  confidence  interval  for  S, 

( - ~~i  >  ®)  -  (15  A,  00) 


*??**  *****$P 


'WtGZs  warns  ■* 


£  I 


<t 

i> 


<<• 

is*. 
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and  the  two-sided  95#  confidence  interval  for  8, 

- i__  ,  . — L_-  j  =  (14.6,  26.5). 

log  -335-  log 


The  best  estimate  for  0  is 

A  i 


log 


a  19.5  months. 


•95 


To  solve  (iii)  we  multiply  the  answers  in  (li)  by  log  2,  Eence 
the  maximum  likelihood  estimate  for  x  K  is  (l9-5){ -693)  ~  13*5  months. 

One  and  two-sided  95S^  confidence  intervals  for  x  ^  are  given  by  (10.7,  od  ) 
and  (10. 1,  18.4)  respectively.  One  can  interpret  (10=7,  00)  as  a  one¬ 
sided  tolerance  interval  in  the  following  sense:  Based  on  the  data  and 
assuming  the  exponential  distribution  we  can  assert  with  95#  confidence 
that  at  least  50$  of  the  items  survive  10. 7  months. 

It  is  interesting  to  raise  the  question:  Suppose  one  knew  the 
actual  failure  times  of  the  50  tubes  which  fail.  How  much  would  our 
estimates  and  confidence  intervals  change?  A  reasonable  assumption  is 
that  the  total  life  of  the  failed  items  is  about.  25  months.  This  amounts 
rough]y  to  assuming  that  the  50  failures  are  uniformly  distributed  over 
one  month.  Thus  T(t#)  *  25  +  950  =  975*  As  a  good  estimate  of  S  with 


very 


little  bias  we  take  *9  *  *  ^15.  „  15,1  months.  One  and  two- 


51 


elded  95#  confidence  intervals  for  6  are  given  by  substituting  in  (17) 
and  (18).  Thus  the  one-sided  95#  confidence  interval  for  9  is  given  by 

(15  A,  00) 


g(HL_ ,  oW-m- ,  m) 

Lfo5<xos>  /  '  / 


and  the  two-sided  95#  confidence  interval  for  0  is  given  by 


w. 

?>.  . 


*  jaaatte.# 
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2T(t») 


X.0£5^102^  X.975^100^ 


(  1950  1950  A 

(  I31TH  »  W&  J 


(14.8,  26.3). 


The  beat  estimate  for  x  K  1b  (19-5)(.6S3)  *  13*5  months.  One  and 

o 

two-sided  95$  confidence  Intervals  for  x  K  ere  (10.?,  a> )  and  (10. 3,  18.2) 

•  y 

respectively.  It  would  appeir  that  little  is  gained  from  actual  knowledge 
of  the  failure  times.  More  will  be  said  about  this  later. 

Problem  7:  20  items  are  tested  one  at  a  time.  If  the  item  falls 

before  1000  time  units  have  elapsed,  the  experiment  is  stopped.  If  the 
item  is  still  living  after  1000  time  units  have  elapsed,  the  experiment 
is  also  stopped,  5  items  are  observed  to  fail  with  failure  times  100, 

400,  600,  800,  900  and  15  items  are  still  living  at  1000  hours.  Give  95$ 
one-sided  confidence  intervals  for  Fr(T  >  t*  •  1000),  the  probability 
of  surviving  t*  *=  1000  tine  units .  . 

Solution  1:  In  the  notation  of  this  section,  n  »  20,  t*  *  1000, 
fj.  m  .05.  A  non-parametric  solution  is  given  by  substituting  in  formula 
(3).  Thus  ve  are  95$  confident  of  the  validity  of  the  assertion  that  the 
probability  of  surviving  t*  *>  1000  time  units  is 


1  +  §  F.05(12>  3°>  1  +  1  (2.09) 


9^3  " 


Put  in  reliability  language,  ve  are  95$  confident  of  the  correctness  of 
the  assertion  that  the  reliability  is  .544  over  the  time  interval 
t*  ■  1000  units. 

Solution  2:  Another  solution  is  obtained  by  assuming  that  the  under* 
lying  distribution  is  exponential  and  applying  (20) .  We  first  calculate 
.  Tt  +  15t*  »  100  +  400  +  600  +  800  +  900  +  15(1000)  *  17800. 
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Bubstitutlng  In  formula  (20)  we  can  be  9556  confident  of  the  correctness 
of  the  assertion  that  the  probability  of  surviving  tine  t#  *  1000  ie 

V  ' 

w  '  %*  X>05(l2)/2T(t*  »  1000)  _  1000  (2i.026)/35600 

d.  c  *  e  »  .>5* 

Put  in  reliability  language,  we  are  9556  confident  of  the  correctness  of 
the  assertion  that  the  reliability  is  '>  -55^  over  the  time  interval 
t*  >  1000  tine  units. 

It  should  be  noted  how  close  the  two  results  (non-parametric  and 
exponential)  are.  Because  of  its  validity  under  much  more  general  condi¬ 
tions,  one  would  normally  prefer  the  non-parametric  solution  1. 
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SECTION  4 

ESTIMATES  OF  BOUNDED  RELATIVE  ERROR  FOR  0 

Problem:  To  give  an  estimation  procedure  for  the  mean  life  0 
having  a  small  relative  error.  Put  more  precisely,  give  a  procedure 
which  will  yield  an  estimate  which  is,  with  some  preassigned  confidence 
1  -  oL,  within  a  certain  percentage  (100  <f  percent)  of  the  true,  but 
unknown  mean  life  0.  In  practice,  oi  and  cf  will  usually  be  snail. 

Approximate  Solution:  In  the  exponential  case,  the  answer  involves 
finding  r,  the  number  of  failures,  such  that 


Such  a  requirement  will  in  general  make  it  necessary  that  r  be  large. 
Let  T  be  the  total  life  associated  with  observing  r  failures  and  let 

X 

©r  -  Ty/r.  Then  it  can  be  assumed  safely  that  ^r(0r  -  0)/8  is  approxi¬ 
mately  distributed  as  N(0,  1).  Thus  to  meet  the  conditions  imposed  by 
equation  (l),  means  that  r  must  be  chosen  in  such  a  way  that 

(?)  r  i  2 


where 

c*.  -  2.576 

if 

a C  *=  .01 

" 

-  1,960 

oL  “  -05 

-  1,645 

oi.  ■  .10  . 

If  4  m  .01,  .05,  -10  and  cL  »  .01,  .05,  10  the  values  of  r  required 

are  tabulated  in  Table  6. 


•3  AT- 


m 

m 


*#  V 

i 

+e* 

f 
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f 

i 

v 


t 


TABLE  6 
Values  of  r 


.01 

-05 

.10 

.01 

66,4 00 

38,400 

27,100 

.  -05 

2654 

1537 

1082 

.10 

664 

384 

271 

Remark:  The  exact  solution  to  this  problem  involves  considerations 
analogous  to  those  .in  the  paper  "Estimates  of  Bounded  Relative  Error  in 
Particle  Counting"  by  M.A.  G±T3hick,  H.  Rubin,  and  It.  Sitgreaves  in  the 
ANNALS  OF  MATHEMATICAL  STATISTICS  26,  276-285,  1355 •  The  values  of  r 
obtained  in  the  range  0  c  ofc  5  .10,  0  <  <f  —  .10  are  almost  identical 

with  those  tabulated  above.  Further  the  "best"  estimator  of  9  in  a 
minimax  sense  for  fixed  r  corresponding  to  Vue  loss  function, 

(3)  L(@,  a)  =  0  if  3  -  /  <  |  i  1  +  <T 

*  1,  otherwise 

is  given  by  the  estimator 


<*> 


?.£T_ 


a  «= 


r  ioe  Hf 


e  _  a  *r 

However,  for  0  <  a  *  .1.  a  —  9^  -  — -  ,  since 
/  '  r  r 

T  *2  «  v 

»  *  y-  (l  -  ^~r  *  °(/  ))  • 


I 


X. 


RerrarK :  One  nan  show  that  ovr  confidence  in  the  validity  of  the 
assertion  that  1  -  <j  •£.  |  <  1  ■*  <f  where  a  .*  2cf  Tr/r  log  ~~  , 
is  given  for  any  ureas signed ..  r  by 

/"tl  +  ^  log 


(5)  ,  Pr(l  ~  £  <  §  <  1  +£)  «, 


x**V* 

j. — ^ 


j  n  „  <M  JL  i  ,-vt  /  j-  *^£\ 
®  '  2k'  ‘tCto  \1  -■  £  • 


-  |1  P”k;  (i  'i]  4  i'°s(W)J  “  £  p&;  (1  -  ^  $  toefj—fil 


For  example,  choose  r  •  10  and  <4  «  10,  then  it  is  readily  verified 

that 

(1  *  S)  log(~~™-)  -  11.1)  ~  log(y)  <*11 

sind 

(l  --  <H  f  log(.i~t|l)  *>-9. 

Renee  (5)  becomes 

CP  qp 

Prfl  ~  $  <  |  <  l  *  £)  »  J2t  P^k:  11  ^  "  21,  P(*S  S) 

k=10  k~10 

*  .6595  -  .1*126  =  .21*69. 

In  other  words,  we  can  have  approximately  254  confidence  in  our  assertion 
that  a  ~~  T.^/lO  is  vTcbin  10$  of  the  true  naan  life  9. 


numerical  Bxaeples 

1.  Bov  many  tubes  should  be  tested  In  order  that  there  1b  a 
probability’  of  at  least  .90  that  the  estimate  is  vlthin  10*  of  the 
true  mean  life? 

Solution:  In  the  notation  of  (l)  and  (2),  B  »  .1,  a  •  .1,  and 
cQ  ■  1.645.  Therefore  the  number  of  tubas  tested  should  be 

ft 

>  (lOO)(l.645)  «•  271.  If  the  underlying  distribution  is  exponential 
this  mean*  that  vo  must  observe  at  least  271  failures  in  order  to  get. 
an  estlaate  such  that  ve  can  be  90*  confident  that  the  estimator  is 
vlthin  10*  of  the  true  hut  unknown  mean  life. 

2.  We  have  available  information  from  a  life  test  in  which  5 
failures  occurred  with  associated  total  life  Ty  «  1000.  Assuming 

an  exponential  distribution,  find  the  minimax  estimator  a  associated 
with  the  loss  function 

L(Q,a)  •  0  if  .8  <  g  <  1.2 
■  1,  otherwise. 

Also  coopute  the  confidence  that  we  will  have  in  the  correctness  of 

the  assertion  that  .8  <  ^  <  1.2. 

Solution;  Fran  (4),  the  winlur  estimator  of  0  based  on  the 

5  failures  is  given  by  a  a  ■  2(  ,2)(200)/log(ii|)  ■  80/log(l.5)  *  197. 

To  find  the  confidence  in  our  assertion  that  .8  <  a/0  <  1.2,  we  Use 

formula  (5).  ibis  gives  us 

00  op 

Confidence  •  >  p(k;6.08)-  p(  k;  4.06  J-.7255-.  3829*..  3^6. 

*5 


-3.50- 

SECTIOH  5 

THE  WO  PARAMETER  EXFCMBWTIAL  DISTRIBUTION 

It  baa  been  found  in  mny  problems  of  life  testing  that  there  are 
occasions  when  a  two  parameter  exponential  distribution  is  more  appropriate 
for  fitting  life  test  data  than  is  a  one  parameter  distribution.  By  a 
two  parameter  exponential  distribution  we  mean  a  density  function 
f(x;  0,  A)  such  that 

(1)  f  (x;  ©,  A)  =  ~  e"*x  "  A^9,  x  2  A  >:  0,  9  >  0. 

A  «A.n  be  thought  of  as  a  {guarantee  period  within  which  no  failures  con 
occur  or  as  a  minimum  life.  If  A  =  0,  equation  (l)  reduces  to  the  one 
parameter  exponential. 

Problem:  A  sample  of  n  items  is  drawn  at  random  from  a  population 

whose  p.d.f.  is  described  by  (l).  The  experiment  is  terminated  as  soon  as 

the  first  r  failure  times  x,  &  x0  <-  . . .  <£  x_  become  available. 

jl  c  —  r 

Items  which  fail  are  not  replaced .  Give  "beBt,"  estimates  for  the  unknown 

parameters  A  and  6. 

Solution:  It  can  be  shown  that  x^,  the  time  to  observe  the  first 
failure,  and  T(xr  -  x^),  the  total  life  observed  in  the  interval  (x^ 
are  mutually  independent  and  Jointly  sufficient  for  estimating  A  and  9. 
Sufficiency  means  roughly  that  x^  and  T(xJ,  -  x^)  Jointly  contain  all 
of  the  relevant  information  for  estimating  A  and  0  that  can  be  obtained 
from  the  first  r  failure  times,  xi  -  x2  *  •••  —  XT’  0681  ostimates 
for  A  and  9  in  the  sense  that  they  are  unbiased  and  minimum  variance 
are  given  by 


\ 

i 


S**** Sfcflft- '4»t*  t,, .  \ 


V 

«> 


■*3  !>i 

(2) 

and 

(3) 

where 

(4)  T(xr  -  Xj)  »  in  ~  l){Xg  -  x^)  +  (n  g)(x.j  -  x^j  +  .  .. 

+  (n  -  r  *  X)(xr  -  xr-1V  >  ,  ■  '  '" 

=  ~(n  -  l}x?  *  x£  +  .  -  +  xr>_1  -  (a  -  r  +•  l)*^ 

It  is  often  convenient  in  (?,)  end  (V  to  use  the  fact  that 
?(>r  -  Xx)  a  r<.'(x?)  -  -  1?{xr)  -  nx^, 

where 

TUr)  ”  12  \  +  (n  - 

Confidence  limits  for  Q  are  easy  to  obtain  from  the  fact  that 
2(r  «  1  js/ 9  2T.’(xr  »•  z^}/©  is  distributed  as  "jCfsti'  ~  2^  Thus  for 

r  2s  2,  one  ana  two-sided  100 fl  -  dL)  percent  confidence  intervals  for 


A  -  }•  -  — 

1  n 


0  --  T(xr  -  xxVr  -  1, 


8  are  given  respectively  by 


(5)  ( 

/ 

r  2{r  ~  1)9 

"  2) 

and 

(6)  / 

/  /v 

JlLile 

9 

l-ta-  •  a; 

V  ■  2 

or 

(  zi(*T  -  ^ 

(  *  2 

\  ?  .  : 

2T(x 


V  \ 
*1' 


,xi(^  -  g) 


XI 


*  (2r 
2 


2T(xr 

X?Jc- 

x  2 


*i) 


>  V 

2}  / 

/  , 


*fr  .«4TVUMW<M9 


t 
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it  follows  that  the  desired  100  /*  percent  confidence  interval  for  A 
is  given  by 


(10) 


zr 


3jr  -  l) 

n 


*<**•  -  *!> 

2j!  .  - - 


Since  Zy  b  Wy/r  -  1  we  have,  of  course,  the  same  confidence  interval 
os  before.  However,  is  computable  for  any  r  and  any  V. 

Remark  2: 


Z 


nOc^  -  A)  n^  -  A) 
(r  -  1)6  =  T(xr  -  xx) 


can  be  Interpreted  as  the  ratio  between  the  total  life  between  time  A 
and  x^,  the  time  when  the  first  failure  occurs,  and  the  total  life  be¬ 
tween  x1  and  xr  Inclusive.  Clearly  one  wants  to  reject  the  hypothesis 
that  A  o  o,  if  nx^/T(xr  -  x^)  is  too  large.  It  should  be  noted  that 
under  the  hypothesis  that  A  =  0,  nx^T(?:r  -  X-^)/r  ~  H  ~  2r-2). 

Remark  3:  Either  formula  (7)  or  its  equivalent,  formula  (10),  can 

be  used  to  test  whether  or  not  A  differs  significantly  from  zero.  If  < 
$  T(x  -  x.) 

Xi-wg-  or  equivalently  ~ -  are  >  0,  then  A 

Is  significantly  greater  than  zero  at  the  (l  -  y*)  level. 

Remark  4:  The  100  percent  confidence  interval  for  A  can  be 

interpreted  as  a  one-sided  tolerance  interval.  More  precisely  we  can 

make  the  statement  that  all  items  live  longer  than  x,  -  zr,‘9  ^  •“ 

.  «*-  *  *i>  1  r 

(or  x^  -  zf . -  ■  -  )  with  confidence  V .  100  y* percent  of  these 

assertions  will  be  correct. 
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!  Numerical  Bxamples 

1.  20  items  are  placed  on  test.  Testing  is  terminated  after  one 

has  observed  the  first  10  failures.  Suppose  that  the  first  failure  occurs 
520  hours  after  the  experiment  starts.  The  total  life  observed  between 
the  time  when  the  first  failure  occurs  and  the  time  when  the  tenth  failure 
occurs  is  12000  item  hours.  Assuming  that  the  underlying  distribution  is 
exponential  do  the  following: 

(i)  Test  vrhether  A  >  0  at  the  .05  level. 

(ii)  If  A  >> 0,  find  the  shortest  95$  confidence  Interval  for 
A  and  an  unbiased  estimate  for  A. 

(iii)  Find  an  unbiased  estimate  for  ©  and  one  add  two-sided 
confidence  intervals  for  ©. 

Solution:  (i)  Suppose  that  A  «.  0,  then 

distributed  os  F(2,  18).  From  the  data 


nx. 


*<*io  - 


is 


_^± _ _  K  (ffl)qgp) 

T(xi0  -  xj/9  .  12000/9 


7.80  . 


But  the  upper  5$  point  for  F(2,  18)  i6  3*55  •  Hsnce  A  is  significantly 
different  from  zero  on  the  .05  level.  As  a  matter  of  fact,  since  the  upper 
.5$  point  for  F(2,  18)  is  7*21  and  the  upper  .1$  point  for  F(2,  18)  is 
IO.39,  A  is  significantly  different  from  zero  at  between  the  .001  and 
.005  levels. 

(ii)  From  the  data  -  T{xlQ  -  x^)/9  «  12000/9  ■  1333*  Hence  an 
unbiased  estimate  for  A  is  given  by  x^  -  ©/n  =  520  -  1333/20  a  520  -  67.7 
a  452.3. 


1.  B.  Epstein  and  M.  Sobel,  "Some  Theorems  Relevant  to  Life  Testing 
from  an  Exponential  Distribution,"  Annals  of  Mathematical  Statistics 
2£,  373-381,  1954. 

2.  B.  Epstein,  "Simple  Estimators  of  the  Parameters  of  Exponential  Dis¬ 
tributions  when  Samples  are  Censored,"  Annals  of  the  Institute  of 
Statistical  Mathematics  8,  15-25,  1956. 


f, 


» 


Appendix  3A 


The  material  in  Section  1  cf  Chapter  3  dealing  vith  point  and 

Interval  estimates  for  the  mean  life  9  is  given  in  detail  in: 

3.  Epstein  and  M.  Scbel,  "fome  tests  based  on  the  first  r 
ordered  observation*  drawn  from  an  exponential  distribution," 
Stanford  University  Technical  Report  No.  6,  Wayne  University 
Technical  Report  No-  1,  March  1952 

and 

I'..  Epstein  and  M.  Sobel,  "life  Testing,"  Journal  of  the 
American  Statistical  Association  48,  486-502,  1953. 


Appendix  3B 

Proof  of  the  Theorems  in  Chapter  3.  Section  1 

.111  order  to  show  th*:t  s',  as  given  by  (7)  is  the  maximum  like- 

r,u 

lihood  estimate  we  write  down  the  p.d.f.  of  the  first  r  out  of  n 
ordered  observations  x^a,  ^  ,  ....  *r_n  .  This  is  given  by: 


f{2L  -  j  •  •  •  •  X  )  sS  y-»— 0 

^  l,n'  2,n’  >  r,n;  (n  “  fJT  Qr 


-[£  V  *  <“  -  r>«r,n> 


n  1 

Tr^r)T  0r 


T  /9 


06  X.  x_  „  <a  .  . . .  ■«:  x _ <  co 


in  the  non -replacement  case  and 

„  _  -nx  /9  -  T  /8 

. V»>  >  <’/  *  1  •  <5>r  '  ^ 


t 


'****' *4, 
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in  the  replacement  case.  It  is  very  easy  to  verify  that  the  maximum  of 
f  occurs  at  ■$* »  Ty/r  and  this  proves  (1) . 

The  sufficiency  of  the  estimate  can  be  verified  at  once  by  using 
a  result  in  Crairifr  j[p.  488j  since  f(x^Q,  *2,n*  Xr#n'  can 
factored  as 

*2,n‘  XT,n'  8)  *  *>  hlxl,n'-xi,n’  "•>  *r,n> 


where 


h(x,  x^  x„  _)■  1  0  jx,  „  f  ...  jx 

CL  jLx  1  |U 


oo 


0  otherwise. 
A 


We  next  show  that  the  p.d.f.  of  9  is  given  by  (6).  To  do 

r,n 

this  ve  introduce  r  new  random  variables  defined  as 

yx  -  rnc1  and  -  (n  -  1  +  l)  ^  -  x^),  2  ±  i  £  r 

in  the  non-replaceme nt  case  and 

yx  -  nxx  and  y±  -  n^  -  x^),  2  <  i  £  r 

in  the  replacement  case.  We  can  now  state  the  following  lemma. 

Lemaa :  The  random  variables  y.  defined  above  are  mutually 

■  i*n 

independent  with  ccamon  p.d.f.  ^  e  1  x  >  0  . 

Proof:  In  both  the  replacement  and  non-replacement  case  the  Joint 

p.d.f.  f(x1#n,  x2#a,  ...,  xr#n)  becomes 

r  _ 

-  Zj  jJ9  r  -  y  /© 

«<yl,n>  y2,n*  •  yr,n>  -  8  *  JJ#  1  ,  yt  -c  oo, 

i  *■  1;  2j  •  •  •  #  3P 

tnd  clearly  the  leone  i»  established. 
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Rewriting  ©_  _  in  tern  of  tbs  y.  ,  (l)  become* 

JT  pi*  XpU 

*r,n  •  IVr- 
I*! 

1  •v/fi 

Since  the  characteristic  function  of  the  p.d.f.  -g  e  '  $  x  >  0  is  given 


by  0v(t)  >  (l  -  ite)"1  it  follows  at  once  from  tbs  independence  proved 


in  the  lemma  that 


0g(t)  *  7T0v  (t/r)  -  (1  -  ^rr. 
*  i-l  yi 


From  the  uniqueness  theorem  for  characteristic  functions  one  gets  by 
inversion  that  the  p.d.f.  of  ©r  ^  is  given  by  (6) } 

<6>  fr«  -  F^-TTT  <i>r  e’rjr/®  <  y  »  0 

*  0,  elsewhere . 

To  complete  the  proof  of  the  theorem  in  Section  1  we  show  that 
a 

e  is  unbiased,  efficient,  and  minimum  variance. 
r,n  ' 

Unbiased neas  is  immediate,  since 


E(e) 


r 

-  E(  yjr)  -  re/r  .  e. 


For  efficiency  and  minimum  variance  let  us  c depute  the  Cramer -Rao 
lower  bound 

1  c  V* 

SEE8.’  “7  * 


_  |  p 

with  C  ■  1x1  uon-replacemont  case  and  C  -  n  in  tbs  re¬ 

placement  case.  '-.Thus  log  f  ■  log  C  -  r  log  0  -  Tj/9* 
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Thus 


^  log  f  -  -r/9  +  Ty/©2 


T-  2 


2  2rT  T 


e 


e3 


%  [l  -  2  +  *^{©2/r  +  e^l  .  r/e£ 


0  e 

2 


Hence  the  Cramer -Rao  lover  bound  Is  © /r. 


But  Var  ©_ 


Var  y 


©  /r  and  since  the  assumptions  needed 


for  the  derivation  of  Cramer-Eao  lower  bound  are  clearly  met  in  the  present 

A 

problem,  ©y  Q  is  minimum  variance  and  efficient  since  any  other  estimate 

2 1 

has  variance  at  least  equal  to  ©  /r.  Thus  the  theorem  in  Section  1  is 
completely  established » 

Remark:  It  is  of  interest  to  note  that  while  6  is  "best”  among 
all  unbiased  estimators,  it  is  not  "best"  or  "admissible"  if  one  uses 
other  criteria.  Using  the  language  of  decision  theory,  let  us  consider 
the  loss  function 


L(©>  a)  »  (9  T..*)  , 


where  9  is  the  true  but  unknown  value  ve  are  trying  to  estimate  and 
where  a  is  our  estimate  of  ©  based  on  knowing  the  first  r  failure 
times.  He  would  like  to  choose  the  estimate  a  in  such  a  way  that 

fc(e,  *)]  is  made  as  small  as  possible  in  the  minima*  sense.  It  can 


be  verified  readily  from  results  in  Chapter  11  of  Blackwell  and  Olrshick's 
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book,  "Theory  of  Games  and  Statistical  Decisions,"  tbat  the  beat  choice 
for  a  la  given  by 

~  rS  *<S> 

» -  *  -  ffx  ■  m  • 

/V 

If  the  estimate  a  «  0  is  used  then 

V<8- s>  i  rh.  * 

.  A.  . 

whereas  if  0  i s  used  as  the  estimate  of  0,  then 

V(®>  ®)«F- 

Hence 

£^(0,  "0)/^t(O,  ^)  «  JTTX 

and  one  always  gets  a  waller  expected  loss  by  using  the  estimate  0 
rather  than  3.  Stated  in  the  language  of  decision  theory,  “©  is  an 
"admissible  ainimax"  estimate  far  the  above  loss  function,  while  3  is 
not  admissible.  Here  is  a  case  where  one  does  better  using  the  "biased” 
estimate  *0  rather  than  the  "unbiased"  estimate  *3. 


Appendix  3C 

We  have  seen  that  a  100(1  -  •L)  percent  on  sided  confldenoe  inter¬ 
val  for  the  quantile  x^,  where  Xp  is  the  solution  to  Pr(X  ar  x^)  -  p 
(i.e.,  x^  «  0  log  j)  is  given  by 

2r8r.n  I 
— ;  A~  ■  iiJ  ,  od 

*>■> 


A 


V 

i 

a 

■% 

t 


i- 


* 


S- 

X- 


» 


+■ 


and  that  this  implies  the  tolerance  statement  that  ve  oan  be  100(1  -Ot) 
percent  confident  of  the  assertion  that  the  fraction  of  itene  surviving 


r 


&•  9  iogg 


it  2p. 


The  proof  of  this  assertion  Is  now  given.  Ve  can  be  100(1  -«L)  percent 
confident  of  the  assertion  that  (T  £  x^  <.  ®  ) .  But 
Pr(X  2  T)  »  Pr(X  Xp)  -  p.  Combining  the  last  two  statements  we  can 
say  that  we  are  100(1  -  d>)  percent  confident  that  the  fraction  of  items 
surviving  T  time  units  is  2  p.  And  this  is  what  we  wanted  to  prove. 


Appendix  3D 

It  is  interesting  to  coopers  the  material  in  Section  1  end  Section  2 
of  Chapter  3  in  the  replacement  case.  We  assume  that  one  starts  the  life 
test  at  time  tQ  ■  0  with  n  items  and  replaces  failed  items  at  once  by 
new  items.  In  the  situation  treated  in  Section  1,  the  life  test  is  con¬ 
tinued  until  a  prescribed  number,  r,  of  failures  have  occurred,  and  one 
stops  testing  at  the  random  time  x^Q  (measured  from  the  beginning  of 
time).  The  total  life  observed  up  to  and  including  XT  &  is  the  random 
variable  Tr  ■  nxr  Q  .  In  Section  2,  the  life  test  is  terminated  at  a 
preassigned  time  t*  and  the  number  of  failures  r  that  occur  is  a 
random  variable.  The  total  life  observed  is  prsasslgned  and  given  by 
T*  -  at*.  To  sum  up;  in  Section  1,  the  number  of  failures  is  fixed  in  ad¬ 
vance  and  it  is  the  waiting  time  (and  hence  total  life)  until  the  r'th 
failure  which  is  random;  in  Section  2,  the  time  (and  hence  total  life) 
allotted  to  the  life  test  is  fixed  in  advance  and  it  is  the  maker  of 
observed  failures  that  is  random. 


1-3 
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kl 


r 

-*» 


Point  estimation  of  ©  in  the  case  where  the  number  of  failures, 
r,  is  fixed  In  advance,  Is  very  slx^le.  The  estimator  0r  n  -  nxr  Q/r 
«  Ty/r  is  among  other  things  a  maximum  likelihood  and  unbiased  estimator 
of  0.  However  in  the  case  where  t*  Is  fixed  as  in  Section  2, 
nt*/r  ■  T*/r  is  a  nmximua  likelihood  estimator  of  0,  but  it  is  biased 
(and  in  fact  meaningless  when  r  »  0).  As  a  natter  of  fact  it  can  be 
proved  that  no  unbiased  estimator  of  0  exists  in  this  case.  If  we 
know  apriori  that  nt*  «  T*  >  >  0,  then  an  almost  unbiased  estimator  for 
0  is  given  by 


0* 


nt*  T» 
r  +  1  "  r  +  1 


This  arises  fron  the  fact  that 


-B 

B(0*)  -  0[l  -  e  J, 


since 


[.**/•  - 1] 

In  any  case,  one  can  find  a  point  estimate  of  0  by  solving  the  equation 
0*  «  0  [l  -  e”  numerically. 

When  one  is  dealing  with  confidence  interval  estimation,  tbs 
situations  in  Sections  1  and  2  compare  as  follows: 


) 
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Fixed  r,  random  Tr  Fixed  T»,  random  r 

100(1  -ot)  percent 
confidence  interval. 

One  sided 


Two  sided 


It  is  very  interesting  to  see  that  there  is  a  striking  similarity  even 

though  the  two  situations  are  radically  different.  It  is  curious  that 

2 

only  the  degrees  of  freedom  for  )(.  need  to  be  changed  as  indicated  above 
when  one  goes  from  the  situation  in  Section  1  to  the  situation  in  Section  2. 

Remark:  It  is  interesting  to  note  that  in  the  case  where  T»  is 
fixed,  Cox  (l]  has  given 

2 T»  a f» 

7*< *  *  l> '  W-  f  (2r  - 1) 

as  on  approximate  two  sided  100(l  -  ot)  percent  confidence  interval. 


(  f r  „  ) 

f  2P* 

) 

/  « r  \ 

f  2T*  2T»  ^ 

tT 


Appendix  3B 

We  should  like  to  verify  that  equations  (l6)  snd  (17)  In  Chapter  3* 
Section  2,  Remark  7  generate  100(1  -  oL)  percent  one  sided  confidence 
Intervals  when  data  arise  from  a  truncated  replacement  procedure 
sln(T  jO.  let  us  first  consider  the  case  where  rQ  ■  1.  In  this 

va 

case  if  no  failures  occur  by  time  t*,  we  stop  life  testing  snd  accord¬ 
ing  to  (l6)  give  i  °°)  u  the  one-sided  100(1  -  ct)  percent  ooo- 

fidence  interval.  If,  however,  a  failure  occurs  at  time  d.  t*,  we , 
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\  ,  ' 


•top  the  Ilf*  tut  and  according  to  (17)  give  (  ■  g  ■  ■  ,  oo )  u  the 
one -aided  100(1  -el)  percent  confidence  Interval.  We  wish  to  verify 
that  this  la  true.  This  means  that  vre  wish  to  prove  that  our  assertion 
that  6  is  contained  in  the  system  of  confidence  Intervals  (l6)  and  (17) 
la  correct  with  probability  2  1  •  no  matter  vhat  0  Is.  This  Is 
particularly  easy  to  do  for  the  case  rQ  »  1.  In  this  case  one  can 
summarize  the  results  In  the  following  table. 


Value  of _ 9 


ronl 


Probability  that  confidence  statements 
based  on  (16)  and  (17)  are  correct 


xi«) 


l-o*. 


If  8  >  . ,  then  no  matter  vhat  happens  our  assertion  is  correct. 


If  9^  — 5 - ,  then  our  confidence  interval  will  not  include  0  if 

.x>) 

the  failure  occurs  after  time  t  ■  — t*.  But 


-X®(2)/2 


ProbCr^  >  1 1  9)  m  e  >  <U  Bence  the  probability  that  our  con¬ 

fidence  Interval  does  not  Include  0  is  equal  to  oLt  and  the  probability 
that  our  confidence  statement  is  correct  la  equal  to  1  -  at,  it  rQ  -  2, 
(16)  and  (17)  give  the  following:  If  no  failures  occur  in  (0,  t*),  stop 
the  life  test  and  give  (  ■2gt*  ,  oo)  aa  the  one  sided  100(l  -<*)  percent 

confidence  interval;  if  only  one  failure  occurs  in  (0,  t*)  stop  the  Ilf* 


test  and  give  (~  ~  co)  as  the  one  sided  100(1  -ot)  percent  confidence 
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interval;  if  2  failures  occur  at  time  7,  ^  t*  then  the  appropriate 

£nT0 

100(1  -ct)  percent  confidence  interval  is  given  by  (r—g — -  ,00).  Again 

we  wish  to  prove  that  our  systen  of  confidence  Intervals  is  correct  with 
probability  i  1  -  ^  no  natter  what  9  is.  It  can  be  verified  that  this 
is  so  and  the  results  can  best  be  sumnarized  in  the  following  table. 


3. 


Value  of  0 
2nt* 


0  > 


*>> 

2ttt» 


2nt* 

x5« 


C  9  £ 


2nt* 

?£<£> 


r  »  2 
o 


Probability  that  confidence  statements 
based  on  (16)  and  (17)  are  correct 


e 

c2 

c*«- 


1  -  at- 
•nt*/  9 


where 


-X2w/2 


-at*-/ 9 


£  06 


9  £ 


2nt» 


1  -  oL 


If  rQ  »  3,  one  gets 


Value  of  9 


Probability  that  confidence  is  fcawnents 
baaed  on  (l6)  ar.d  Qt)  are  correct 


9  > 


2nt» 

%> 


a  _  2nt* 
9  -  — 27~ 


1  ~  oL 


2nt*  «  _  2nt* 

XI(2> 


^  e  <4 


-  V  2 


X>) 


*■  1 

1.1  * 


w  «.  — ^ - 

thv 

oL 

For  general  rQj  one  gets 


Value  of  6 


X>) 


X5*> 


T#2*  *  2> 


<.  9  ■£.  2nt» 

x5« 


*>o> 


*M*?#***WS*r  ,.*^.vn«.w#.v  -rr1-,.', 
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In  Chapter  III,  Section  3*  we  dealt  with  life  teat  situation* 
in  which  n  items  are  placed  on  teat,  where  testing  la  discontinued 
after  a  fixed  time  t*  has  elapsed,  and  where  items  which  fail  are  not 
replaced.  In  the  first  part  of  this  section  we  gave  estimation  procedures 

i 

which  depended  only  on  r,  the  observed  rnnber  of  items  failing  in  [o,  t*J, 
and  not  on  ^  f  ^  *  ...  £fr  £  t*,  the  actual  failure  times.  More 
precisely,  we  gave  non-parametric  one  and  two  sided  confidence  Intervals 
for  the  probability  of  surviving  for  a  length  of  time  t*,  and  in  the 
special  case  where  the  underlying  distribution  is  exponential  we  were  able 
to  translate  these  Intervals  into  confidence  statements  about  the  naan 
life  e.  — 


Since  the  sufficient  statistic  for  estimating  0  in  this  problem 


given  by  the  pair  (r,  T(t*))»  5^  -f(n  -  r)t*) 


we  know  that  we 


can  make  better  estimates  and  better  confidence  statements  about  0,  if 
we  use  not  only  r  but  also  T(t*) .  To  carry  this  out  in  practice,  how¬ 
ever,  is  not  easy  since  the  c.d.f.  of  T(t*)  is  expressible  only  in  a 
series  of  many  terms.  The  c.d.f.  Is  given  in  S.  Taksda  and  S.  Shlosda, 
"Statistical  Analysis  of  Life  of  Vacuum  Tubes,"  Hitachi  Review,  pp.  143- 
154,  July,  1954.  See  particularly  page  153*  The  c.d.f.  la  given  by 
equation  2.6.16  in  the  paper.  Thus  one  is  virtually  forced  to  use  approxi¬ 
mate  confidence  intervals  and  to  use  a  certain  amount  of  heuristic  reason- 


One  approach  to  the  problem  of  finding  approximate  confidence 
intervals  is  given  in  the  paper  by  Takada  and  Shlmada  to  which  we  have 


9 
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just  referred.  Essentially,  their  idea  is  as  follows.  If  we  pistes  a 
items  on  test  and  truncate  life  testing  at  tine  t*,  we  can  treat 
T(t*)  m  T.  +  (r  -  r)t*  as  the  sum  3^)  X,,  where  the  X.  are 

53  i-1 

identically  and  independently  distributed  random  variables,  each  possess* 
ing  the  c.d.f. 

P(t)  -  1  -  e’*/®  ,  t  t* 

-  1  ,  t  >  t*. 

Tahada  and  Shiraada  apply  the  central  limit  theorem  to  obtain  an  approxi¬ 
mation  to  the  c.d.f.  of  T(t*)  by  the  normal  distribution.  From  this 
approximation  they  obtain  appropriate  confidence  limits  for  9.  They 
claim  that  this  approximation  is  a  very  good  one  and  give  a  table  which 
states,  far  exsnple,  that  if  ri  *  20,  and  t*/®  ■•Of,  then  an  error 
of  5Jt  is  made.  They  further  state  that  if  n  2r  30,  and  t*/0  2r  .1 
or  n  >  50  and  t*/9  fc  .0$  then  the  arror  associated  with  the- approxi¬ 
mation  Is  lees  than  1$. 

We  have  given  another  approximation  in  equations  ( IT )  and  (18) 
of  Chapter  3.  These  formulae  are  certainly  excellent  approximations  for 
n  large  and  even  for  email  a,  they  should  be  quite  good.  There  ere 
a  number  of  reasons  why  we  believe  that  this  statement  is  correct.  Among 
these  are: 

(i)  If  t*/9  is  small,  then  the  number  of  failures  will  be  small 
and  ths  non-replacement  case  becomes  virtually  s  replacement  case.  One 
can  then  act  as  if  we  were  observing  a  fblsson  process  with  rats  A  » 
for  a  length  of  time  T(t*); 


H|® 


ur®****#  nil 
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<*~7s  - 
A 


A-  2g  jag.*) 


But  It  is  readily  verified  that 


Hence 


O  Iob  f  .  — / »  2 

^r“"r/A  • 


^  '  £<^)  ‘  .a  -  =- M*) 

As  another  way  of  estimating;  A,  we  note  that 
Pr(T  *  t*)  -  1  -  e"  M*  -  £  (£). 

i  •  it# 

Hence  the  statistic  r/n  is  an  unbiased  estimate  of  1  -  e 
and  thus  an  estimate  of  A  is  given  by 

A  -  los 

S\  >v 

As  was  the  case  with  A,  A  Is  also  biased  for  finite  n.  However  as 
n  — *  oo  A  A  also.  Let  us  now  compute  the  asymptotic  variance  of 

O' 

A  •  It  can  be  shown  that  as  n  — *  <l 


Yar  n 

(n  -  r)"]2  (t*)S 


d.'***  <1 


.»«*  a . ***) 


n - -  ---  .  «f  .wse-rs  •/  . . .  smm  '  Tu  r~' 


'JWVWV**-  •'■'•  ■ 


'jffp*.--*-  <-T'-  **»■>  •** V»  »H"- 
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Let  ue  now  compute*  the  ratio  cr^  /<r^. .  it  is  easy  to  verify  that 
-  _ 2/  _2  (At*)2e“At* 

■V^r-  (1:.-*ss  • 


Expanding  the  mmnrator,  ve  let 


(At*)2  e"  Xt*  .  (At*)2  -  (At*)3  +  .... 


And  expanding  the  denominator  ve  get 

(1  -  e~  At#)2  -  (At*)2  [l  -  At*  +  (At*)2^  -  ...1 

g  2 

neglecting  higher  order  terms,  /<r^  becomes 

1  -  At*  ♦  -  . . . 

l  -  At*  +  (A  t*)  j|  •  ... 

It  Is  interesting  that  the  ratio  is  close  to  one,  (i.e.,  A  is 
almost  as  efficient  as  ^  )  particularly  if  At*  is  £  jfi  Indeed,  if 
At*  ■  it  is  readily  verified  that  <r^2/«r^?  ~  •  Although  what 

ve  have  Just  done  is  for  point  estimates,  clearly  similar  results  will  bold 

for  confidence  Intervals.  Also,  it  is  trivially  noted  that  although  ve  vere 

%  \ 
discussing  estimation  of  A,  tbs  conclusions  obviously  apply  to  the  \ 

parameter  &  •  l/A  ae  veil.  The  upshot  of  ths  preceding  discussion 
is  that,  in  case  the  underlying  distribution  is  exponential,  then  the  con¬ 
fidence  intervals  (5)  and  (6)  given  in  section  3,  which  depend  only  on 
the  number  of  failures  r  in  (0,  t*),  are  almost  as  short  as  thorns  based 
on  using  both  r  and  T(t*). 


I 
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It  Is  interesting  to  note  that  if  life  testing  is  terminated  not 
after  a  preassigned  time  t*,  but  after  a  preassigned  total  life  T*, 
then  the  problem  becomes  one  of  malting  appropriate  estimates  of  A  or 
8  a  l/A  when  observing  a  Poisson  process  having  rate  A  -  l/o  far 
a  length  of  time  T*.  Thus  the  considerations  In  Section  2  and  Appendix 
3E  can  be  used,  the  only  difference  being  that  we  replace  at*  in  Sec¬ 
tion  2  by  T*.  We  now  state  a  number  of  results  without  proof. 

Suppose  that  life  testing  stops  after  a  total  life  T*  has  been 

observed.  If  the  underlying  distribution  Is  exponential  with  roan  life  6, 

then  the  number  of  observed  failures,  r,  is  a  Poisson  random  variable 

distributed  vlth  the  probability  law 

_  T* 

Pr(r  .  kje  )  -  p(k;  ~)  -  e  ^  (r*)* k  m  0,  l,  2,  ... 


Using  precisely  the  saaa  arguments  as  In  Section  2,  it  can  be  asserted 
that  if  r  *  number  of  items  which  fail  in  (0,  T*),  than  a  one-sided 
100 (l  -  at)  percent  confidence  interval  for  8  Is  given  by 


- — -  > 

+  2)  J 


and  a  two-sided  100(1-  at)  percent  confidence  Interval  for  8  Is  given  by 


2  T» 

x2 F 

B 


; _  2  T»  | 

*■»'  xi -y*y 


Sots  that  for  r  -  0,  only  one-sided  confidence  limits  sake  sense 


■/t 


ffTiTri^'i —  • 


Another  kind  of  situation  is  where  date  become  available  ae  the 
result  of  the  following  rule  of  action:  Reject  if  rQ  failures  occur 
before  total  life  7*  has  been  used  up;  accept -if  fewer  than  rQ  fail¬ 
ures  occur  by  the  time  one  has  observed  a  total  life  of  T*  (it  is  assumed 
that  tq  and  T»  are  preassigned).  Zn  the  event  that  one  rejects,  ex¬ 
perimentation  stops  at  T(tJ.  ),  the  total  life  observed  up  to  and  includ- 

o 

lng  T'  ,  the  r.'th  failure  time.  In  the  event  that  one  accepts  the 
ro 

total  life  observed  will  be  T». 

Using  precisely  the  same  considerations  as  In  Section  2  and  in 
Appendix  3E,  we  can  assert  that  if  the  number  of  observed  failures  in 
(0,  T*)  is  0  4  k  4  r  •  1,  then  a  one-sided  100(l  -  «fc)  percent  con¬ 
fidence  interval  is  given  by  - 


When  r  -  rQ#  i.e.,  if  T(r^  I*  then  the  appropriate  100(1  -  «*•) 

*o 

percent  confidence  Interval  is  given  by 


Similar  results  can  be  conjectured  for  two  sided  100(1  -  oL)  percent 


confidence  intervals.  Tbs  results  are 


if  k  -  0 


t 
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the  desired  one-sided  ^Of>  confidence  interval  is  [103^00 3  (i.e. »  we  can  toe  90^  confident  of  our 
assertion  that  6  >  103*0  *  i 
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TABUS  4(a) 

Values  of  -■■  -v«— ~  — - —  for  a  «  .01 

1  4  (^)Fa(er+2'2n‘2r) 


0 

1 

2 

a 

4 

5 

6 

7 

8 

9 

1 

.0100 

2  . 

.1000 

.0050 

3 

.2158 

.0588 

.0033 

4 

.3162 

.1406 

.0420 

.0025 

5 

.3981 

.2220 

.1056 

.0327 

.0020 

6 

.4640 

-2945 

.1731 

,084? 

.0268 

.0017 

7 

.5181 

<3567 

.2362 

.1422 

.0708 

.0226 

.0014 

a 

,5622 

.4098 

•  2933 

.1981 

.1210 

.0608 

.0197 

.0013 

9 

.5996 

.4561 

-3435 

.2500 

.1709' 

.1052 

.0534 

.0173 

.0011 

10 

.6309 

.4956 

.3883 

.2971 

.2182 

.1503 

.0932 

.0475 

.0155 

.0010 

11 

.8579 

.5302 

.4280 

•  3396 

.2622 

,1938 

.1344 

.0836 

.0428 

.0141 

12 

.6814 

-5607 

.4627 

•3775 

.3025 

.2349 

.1747 

.1215 

•0759 

.0390 

13 

,7016 

.5871 

•4937 

.4195 

<>3390 

.2730 

.2128 

.1569 

.1103 

.0694 

14 

.7196 

.6109 

.5215 

.4435 

-3724 

.3080 

.2488 

.1947 

.1457 

.1018 

1  15 

.7357 

.6323 

.5469 

.4717 

.4029 

.3404 

.2822 

.2288 

.1793 

.1345 

*  16 

.7498 

..65. 10 

•  5693 

.4 969 

.4309 

.3701 

.3134 

.2607 

.2117 

.1663 

17 

.,7627 

.6683 

.5903 

.5201 

.4569 

.3976 

••3423 

.2907 

.2421 

.1970 

18 

.7742 

.,6838 

.6086 

•5419 

,4803 

.4226 

.3691 

.3187 

.2709 

.2261 

19 

.7848 

.6982 

.6257 

,5610 

.5017 

.4459 

.3937 

.3448 

.2979 

.2538 

20 

.7943 

,7111 

;64l7 

.5790 

.5220 

.4682 

.4175 

•3689 

.3234 

.2799 

25 

.8317 

.7624 

•7042 

.6509 

.6017 

-5562 

.5117 

.4697 

.4294 

♦3902 

30 

.8576 

.7985 

.7480 

.7025 

.6596 

.6194 

.5803 

.5430 

,5077 

.4730 

40 

.8913 

.8453 

.8058 

.7698 

.7360 

.7042 

.6732 

.6434 

.6141 

.5862 

50 

.9121 

.8744 

.8422 

.8127 

.7853 

.7583 

•7326 

.7079 

.6838 

♦f596 

75 

.9404 

.9147 

,8927 

.8722 

.3529 

.8342 

.8166 

■  7991 

.7822 ' 

•7655 

100 

.9550 

•9356 

.9187 

.9030 

.8883 

.8742 

.8605 

.8471 

.8332 

.8212 

150 

•  9697 

.9566 

.9452 

.9346 

•  9246 

-9151 

.9056 

.8965 

.8873 

.8788 

200 

•  9772 

.9672 

.9586 

*9506 

.9430 

.9358 

.9286 

.9216 

.9150 

.9081 

300 

.9847 

.9780 

.9722 

.9 668 

.9617 

.9569 

,9520 

.9473 

,9428 

.9382 

400 

.9685 

.9835 

•9791 

•9751 

.,9712 

.9675 

.9639 

.9603 

-9570 

♦9535 

500 

,9908 

.9868 

.9833 

.9800 

•  9770 

.9740 

•9711 

.9683 

,9655 

♦9627  ~ 

£ 
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n  ,0009 

12  .0128  .0008 


13 

,0358  .0118 

.0008 

— 

14 

.0640  .0331 

.0109 

.0007 

15 

.0944  .0594 

.0307 

.0102 

.0007 

£ * 

.1249  .0878 

.0554 

.0287 

.0085 

.0006 

17 

.1553  .1168 

.0821 

.0519 

.0269 

.0076 

.0006 

18 

.1842  .1453 

.1096 

.0772 

.0488 

.0253 

.0069 

.0006 

19 

.2126  .1733 

-1367 

.1033 

.0727 

.0460 

.0237 

.0063 

.0005 

£0 

.2387  .2000 

.1634 

.1292 

.0976 

.0688 

.0436 

.0223 

,0056 

,0005 

25 

.3520  .3165 

.2817 

.£479 

.2155 

.1846 

.1553 

.1276 

.1007 

.0765 

.0541  | 

30 

.4383  .4056 

.3738 

.3436 

.3131 

.2834 

•2555 

.2281 

.2013 

.1757 

.1509 

40 

•5592  .5313 

.5059 

.4802 

.4541 

.4290 

.4052 

.3814 

.3580 

.3344 

.3122 

50 

.6372  *6146 

.5923 

.3700 

•  5489 

•  5274 

.5066 

.4861 

.4657 

.4457 

.4261  1 

| 

75 

-7*94  -7335 

.7173 

.7020 

.6868 

.6711 

.6558 

.6408 

.6258 

.6112 

•  5970 

100 

.8082  .7958 

•7840 

.7720 

.7598 

.7478 

•7362 

.7247 

.7132 

.7019 

.6906 

150 

.8674  .8617 

.8531 

.8448 

,8368 

.8283 

.8201 

.8120 

.8043 

.7967 

*7888 

200 

.9017  .8953 

.8888 

.8825 

.8762 

.8770. 

.8635 

>8575 

.8515 

.8457 

•8397 

300 

.9339  .9295 

.9254 

.9211 

.9165 

.9123 

.9082 

•  9041 

.9001 

.8962 

.8920 

400 

.9502  .9469 

.9437 

.9405 

.9370 

.9338 

.9307 

.9277 

-9246 

•9216 

.9185 

500 

.9601  .9575 

•9549 

.9523 

.9496 

.9470 

.9445 

.9420 

•  9396 

.9372 

.9347 

S 


Values  of 


3.83 

TAMLE  h(a) 


1  +  F{?.J*2,2n~2r) 

u^I  u- 


for  a  «*  .05 


I 

■— WW'M 

0 

i 

2 

3 

4 

5 

6 

7 

8 

9 

i 

1 

,0500 

H 

£ 

.2237 

,0254 

3 

,3686 

•  1353 

.0170 

4 

,4728 

.2488 

.0977 

.0127 

5 

>5495 

.3425 

»l894 

.0765 

.0102 

6 

..6067 

,  41.8.1 

.2714 

*1531 

c.0629 

.0085 

m 

mn 

7 

,6518 

.4792 

.3411 

.2252 

.1288 

.0546 

.0073 

& 

8 

„68?9 

.5295 

.4000 

,2894 

.1928 

.1111 

.0464 

.0064 

I 

0 

.7171 

,5706 

.,4502 

■  3448 

*2513 

.1689 

*0977 

.0411 

.0057 

1 

10 

.,7413 

=6057 

.4932 

;3933 

.3038 

.2226 

.1499 

.0873 

*0397 

.0051 

m 

» 

11 

,7618 

.6353 

.5300 

*4357 

.3500 

,2710 

.1998 

*1389 

.0787 

.0333 

1 

fe 

12 

.7792 

.6611 

.5618 

.4727 

•  3912 

-3156 

.2451 

•l£o8 

.1230 

.0719 

13 

,7941 

.6834 

■  >893 

*5031 

.4276 

•3552 

.2874 

*2239 

.1656 

.1126 

| 

14 

,8074 

.7035 

.6144 

•  5340 

.4598 

.3906 

.3253 

.2640 

*2058 

*1529 

J  ** 

15 

,8188 

..7209 

>6369 

•3597 

.4889 

.4223 

*3596 

.3003 

.2440 

•1899 

1 

*> 

16 

,8294 

.7360 

*6557 

.3835 

*5161 

.4512 

.3915 

*3333 

,2787 

,2266 

i 

» 

17 

.8383 

.7498 

6739 

.6045 

*5394 

•4785 

.4200 

.3644 

.3106 

.2597 

j. 

* 

18 

.8467 

.7623 

6897 

.6229 

-5611 

,5019 

.4459 

.3923 

*3407 

*2913 

jj 

f 

19 

•8543 

<■7733 

.7042 

,64.10 

,5814. 

,5240 

,4705. 

4178 

*3679 

*3205 

20 

.,8610 

,7838 

*7177 

•  6559 

*5993 

.5447 

.4926 

,4422 

3941 

.3470 

25 

,8872 

,8236 

-7692 

..718O 

,6709 

.6250 

*5807 

.5383 

*4959 

•4558 

30 

,9050 

.8512 

.8046 

.7615 

.7200 

.6812 

.6434 

•6059 

,5705 

*5344 

;  i 

40 

«9279 

.8867 

.8510 

,8174 

.7855 

*7553 

.7248 

.696*1 

,6684 

*6394 

i  1 

50 

,9418 

,9087 

.8794 

•  8523 

.8260 

.8013 

*7770 

•753-1 

.7300 

*7071 

1 

75 

,9608 

*9384 

.9184 

*8999 

.8822 

,8648 

.8481 

*8321 

,8156 

.8003 

100 

.9705 

.9534 

.9385 

.9244 

,9107 

.8977 

.8850 

.8727 

,8601 

.8483 

150 

,9802 

.9688 

*9586 

.9491 

.9401 

*9313 

,9227 

.9140 

*9057 

.8974 

200 

9851 

.9765 

.9688 

.9616 

.9548 

*9482 

*9417 

•9352 

.9290 

.9226 

i 

300 

•  ..9901 

,9843 

9791 

.9743 

*9697 

•  9653 

.9609 

.9567 

*952? 

,9482 

400 

.9925 

,9882 

*9843 

•  9807 

*9772 

.9739 

*9706 

.9674 

*9643 

.9611 

i 

500 

,9940 

*9906 

9874 

.9845 

.9818 

*9791 

*9765 

.9739 

*9714 

,9688 

f- 

* 
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Values  of 


TABLE  4(b)  -  cent. 

_ 1 _ 

1  +  Fa(2r+2,2n-2r) 


for  a  »  .0? 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


11 

.0046 

12 

.0304 

.0043 

13 

.0658 

.0281 

.0039 

— 

14 

.1041 

.0611 

.0260 

.0036 

♦  15 

.1418 

.0965 

.0568 

.0242 

.0034 

*16 

.1779 

.1407 

.0903 

.0531 

.0227 

.0032 

' 

17 

.2119 

.1662 

•  1239 

.0846 

.0499 

.0213 

.0030 

18 

.2441 

.1989 

.1563 

.1163 

.0797 

.0492 

.0201 

.0028 

19 

.2738 

.2296 

.1876 

.1473 

.1099 

•0753 

.0444 

.0190 

.0027 

20 

.3021 

.2586 

.2170 

•1773 

.1394 

.1037 

.0712 

.0422 

.0181 

.0026 

25 

.4166 

•3792 

•3413 

.3053 

•  2793 

.2354 

.2021 

.1707 

.1396 

-1099 

.0822 

30 

.3011 

.4666 

.4334 

.4026 

•3695 

.3390 

.3085 

.2790 

.2497 

.2209 

•1931 

40 

.6132 

.5859 

•5595 

•  5339 

.5075 

.4822 

.4572 

.4329 

.4083 

.3846 

•3615 

50 

.6840 

.6627 

.6409 

.6182 

•  5976 

.5763 

•5556 

.5349 

*5143 

.4943 

.4745 

75 

.7842 

•7694 

.7542 

.7383 

.7237 

.7086 

.6938 

.6792 

.66*7 

.6503 

.6360 

100 

.8361 

.8245 

.8127 

.8014 

.7898 

.7784 

.7672 

.7558 

.7443 

•7330 

.7224 

150 

.8893 

.8817 

.8737 

.8659 

•8576 

.8499 

.8424 

.8349 

.8272 

.8196 

.8121 

200 

.9166 

,9106 

.9045 

.8987 

.8925 

.8866 

.8809 

.8751 

.8694 

.8638 

.8580 

. 300 

.9442 

.9400 

•9359 

•9319 

.9280 

•  9239 

.9199 

.9160 

.9121 

.9083 

.9043 

400 

.9580 

.9548 

.9517 

•  9487 

.9458 

•  9427 

•9397 

.9367 

•  9337 

.9309 

.9280 

500 

.9664 

.9638 

.9613 

.9589 

•  9566 

.9541 

•9516 

.9492 

-9469 

.9446 

.9423 

f 
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TABLE  4(c) 


0 

1 

2 

3 

4 

5 

6 

7 

Q 

9 

1 

.1000 

2 

•3165 

.0513 

3 

.4644 

•  1957 

.0345 

4 

•  5626 

•  3205 

.1426 

.0260 

5 

•6313 

.4158 

.2469 

.1124 

.0209 

6 

.6810 

.4892 

•  3331 

.2011 

.0926 

.0174 

7 

•7194 

.5474 

.4039 

.2786 

•  1695 

.0787 

.0149 

8 

.7498 

•  5942 

.4619 

.3444 

.2395 

.1450 

.0686 

.0131 

9 

.7745 

.6319 

.5102 

.4011 

.3012 

.2100 

.1295 

.0608 

.0116 

10 

.7943 

.6627 

.5502 

.4487 

•  3540 

.2671 

.1873 

•  1159 

.0546 

.0105 

11 

.8112 

.6897 

.5848 

.4890 

.4000 

•  3171 

.2405 

.1695 

.1047 

.0495 

12 

•8253 

,7124 

6?  46 

-  .5245 

.4408 

•3611 

.2881 

.2189 

.1541 

•0955 

13 

.8376 

.7326 

.6403 

•5556 

.4762 

.4007 

•3303 

.2639 

.2002 

.1418 

14 

.8485 

.7497 

.6623 

.5826 

•  5076 

•4355 

.3686 

•  3043 

.2427 

.1852 

15 

.0576 

.7642 

.6831 

.6073 

•5366 

•4673 

.4026 

•  3413 

.2820 

.2256 

16 

.8658 

.7780 

.7000 

.6286 

.5607 

,4959 

,4344 

.3744 

.3176 

•2632 

17 

.8731 

.7897 

.7163 

.6481 

.5830 

.5208 

.4623 

.4058 

.3497 

•2974 

18 

.8798 

.8004 

.7303 

.6661 

.6034 

.5435 

.4878 

•4331 

.3804 

.3285 

19 

.8858 

.8101 

.7430 

,6814 

.6224 

.5645 

,5106 

.4587 

.4071 

•3584 

20 

.8913 

.8190 

.7547 

•  6956 

•  6390 

.5841 

•  5319 

.4815 

.4324 

.3846 

25 

.9121 

.8531 

.8004 

.7516 

•  7054 

•  6596 

.6160 

.5734 

•5313 

.4911 

30 

.9262 

.8765 

.8322 

•7906 

.7508 

.7114 

•6743 

.6381 

.6021 

•5673  , 

4o 

.9441 

.9060 

.8723 

.8404 

•8097 

.7794 

•7510 

.7232 

•6953 

.6685 

50 

.9549 

.9245 

.8972 

.8713 

.8465 

.8215 

•7984 

•7757 

•7530 

•7311 

79 

•9697 

.9491 

.9305 

. .9129 

.8960 

.8790 

.862 9 

•8473 

.8317 

.8168 

100 

•  9772 

.9617 

.9477 

.9344 

•  9216 

•  9086 

.8967 

.8848 

.8728 

.8613 

150 

.9847 

.9743 

•  9649 

•  9559 

.9473 

.9386 

•9305 

•  9225 

•9143 

•9065 

200 

.9885 

'.9807 

.9735 

•  9669 

.9604 

.9538 

,9477 

.9416 

•9355 

•9295 

300 

•9924 

.9871 

.9823 

•9779 

•  9735 

.9691 

•  9650 

.9610 

•9568 

•9526 

4oo 

.9943 

•  9903 

.9867 

.9834 

.9801 

•  9768 

•  9737 

.9706 

•9675 

•9645 

500 

.9954 

.9922 

.9893 

•9867 

.9840 

.9814 

.9789 

•9764 

•9739 

•9715 

( 


TABLE  Me)  -  cont. 

V&luee  of  - api - -  for  a  *  .10 

1  ♦  (~~)  Fa(2t42,2n-2r) 

10  U  12  13  14  1?  16  1?  18  19  20 


11 

.0095 

12 

.0452 

.0087 

13 

.0879 

.0417 

.0081 

.14 

.1311 

.0814 

.0386 

>15 

.1719 

.1220 

.0756 

16 

.2102 

.1605 

.1136 

17 

.2460 

.1969 

.1500 

18 

.2789 

.2312 

.1852 

19 

,3092 

.2628 

.2181 

20 

.3381 

.2923 

.2486 

25 

.4513 

.4119 

,3740 

30 

.5325 

.4962 

.4648 

40 

.6409 

.6137 

.5873 

50 

.7085 

.6863 

.6646 

75 

.8012 

.7860 

.7711 

100 

.8494 

.8378 

.6265 

150 

.8984 

.8906 

.8830 

200 

-923* 

.9175 

.9117 

300 

.9488 

.9448 

•9409 

400 

.9614 

.9584 

.9554 

**500 

.9690 

.9666 

.9642 

-0075 

.0360 

.0070 

.0709 

.0337 

.0066 

.1068 

.0667 

.0317 

.0062 

.1413 

.1008 

.0628 

.0299 

.1750 

.1337 

.0951 

.0595 

.2024 

.1660 

.1264 

.0903 

.3311 

.3014 

.2643 

,2302 

.4323 

.4006 

.3689 

.3373 

.5616 

•  5367 

.5112 

.4855 

.6433 

.6224 

.6008 

.5796 

.7566 

.7424 

.7274 

.7127 

.8155 

.8047 

.7932 

.7818 

.8755 

.8683 

.8605 

.8528 

.9061 

.9006 

.8946 

.8888 

.9371 

.9334 

>9295 

.9255 

.9526 

.9498 

.9468 

.9438 

•9619 

-9597 

•9573 

.9549 

.0058 

.0283 

.0055 

.0564 

.0268 

.0053 

.1963 

.1631 

.1310 

.1001 

.3068 

.2768 

.2472 

.2178 

.4603 

.4356 

.4114 

.3877 

.5587 

.5383 

•  Slfll 

.4984 

.6962 

.6839 

.6698 

.6559 

.7706 

•  7595 

.7486 

•7378 

.8451 

.8376 

.8304 

.8226 

.8831 

.8774 

.8717 

.8662 

•9217 

•  9179 

.9141 

.9104 

•9409 

.9380 

.9352 

•  9324 

•9526 

9503 

.9480 

-9457 

>" 
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TABLE  4(d) 

Values  of  - -gv  — ■  for  a  »•  .25 

1  +  <§?>  y2r*2,2n-?;r) 


0 

1 

2 

3 

*  5 

6 

7 

8 

9 

1 

.2500 

2 

.5000 

.1340 

3 

.6303 

•3268 

.0915 

% 

*7067 

.4559 

.2427 

.06S>4 

5 

-7576 

.5464 

.3597 

.1938 

.0559 

6 

.7937 

.6112 

.4469 

.2964 

.1613 

.0469 

7 

.8206 

.6593 

•5133 

•3788 

.2532 

.1382 

.0403 

8 

.8412 

.6972 

■5666 

.4448 

.3292 

.2203 

.1208 

.0354 

9 

.8571 

•7273 

.6087 

.4983 

.3922 

-2915 

.1949 

.1073 

.0315 

10 

.8703 

•7525 

.6446 

.5418 

.4444 

.3511 

.2607 

.1756 

.0965 

.0283 

11 

.8614 

•7728 

.6742 

-5797 

.4895 

.4016 

.3178 

.2358 

.1592 

.0077 

12 

.8909 

•7914 

.6983 

.6114 

.5263 

•4459 

•3662 

.2900 

.2163 

.1456 

13 

.8990 

.8065 

.7208 

.6378 

•  5590 

.4825 

.4093 

•3366 

.2676 

.1990 

14 

.9056 

.8186 

•7394 

.6627 

.5882 

.5172 

.4465 

.3783 

.3120 

.2475 

15 

.9119 

.8304 

.7545 

.6834 

.6128 

.5453 

.4799 

•4153 

.3517 

.2699 

16 

.9170 

.8408 

.7692 

.7019 

.6349 

.5723 

•  5093 

•4477 

.3877 

.3286 

17 

.9216 

.8496 

.7825 

.7172 

.6549 

.5952 

.5354 

•4775 

.4195 

•3636 

18 

.9257 

.8575 

.7937 

.7324 

-6731 

.6161 

.5588 

.5035 

.4489 

.3947 

19 

.9295 

.8647 

.8039 

.7454 

.6897 

.6352 

.5803 

.5271 

.4755 

.4237 

20 

.9328 

.8712 

.8132 

.7573 

.7039 

.6510 

.6000 

.5488 

.4991 

.4508 

25 

.9459 

.896I 

.8490 

.8036 

•7599 

.7179 

.6748 

.6331 

.5923 

•5525 

30 

.9548 

.9130 

•0733 

.8355 

.7985 

.7622 

•7262 

.6908 

.6560 

.6219 

40 

.9659 

.9340 

.9040 

.8755 

.8475 

=8195 

•7924 

•7656 

.7390 

.7128 

5C 

.9726 

•  9469 

.9228 

.8998 

.8771 

.8547 

.8327 

.8110 

.7893 

.7679 

75 

.9817 

-964j 

.9482 

•9326 

•  9173 

.9025 

.8877 

.8730 

.8584 

.8440 

100 

.9862 

-9731 

.9609 

.9491 

-9375 

.9263 

.9151 

.9039 

.8927 

.6818 

15C 

.9908 

.9819 

•9737 

.9658 

.9580 

•9504 

•  9429 

.9353 

.927 8 

.9204 

200 

•9930 

.9864 

.9802 

.9743 

.9664 

•9627 

•9570 

.9513 

.9456 

.9400 

300 

.9954 

.9909 

.9868 

.9828 

•9788 

•9750 

•9712 

.9674 

•9635 

.9596 

400 

•9965 

•  9932 

.9901 

.9871 

.9841 

.9812 

.9783 

.9755 

.9726 

•9697 

500 

.9972 

•9945 

.9921 

.9896 

.9873 

.9850 

.9626 

.9803 

.9780 

.9758 

TABLE-  4(4)  -  «»*• 

Valuae  of  - ~ - 7  tor  a  * 

1  ♦  (~)Fa(2r+2,2n-2r) 


13  14  15  16  17  18 


i 

11 

.0236 

u 

1 

12 

.0801 

.0237 

I 

1 

13 

.1342 

.0742 

.0219 

1 

1 

14 

.1852 

.1250 

.0689 

.0203 

1 

i 

£  15 

.2302 

.1724 

.1165 

.0643 

.0190 

s 

4 

16 

.2706 

.2151 

.1563 

.1091 

.0602 

.0178 

& 

jg. 

17 

.3095 

.2551 

.2019 

.1515 

,1026 

.0567 

.0168 

id 

.3419 

.2912 

.2402 

-1913 

.1429 

.0968 

.0535 

.0158 

'» 

;  19 

'  *3731 

-3241 

-2792 

.2281 

.1817 

•1351 

.0916 

.0307 

.0150 

4v.« 

i*  . 

20 

.4015 

.3538 

.3074 

.2577 

.2162 

.1715 

.1282 

.0870 

.0482 

.0143 

- 

25 

.5129 

.4730 

•4335 

.3949 

•  3571 

.3181 

.2602 

.2449 

.2072 

.1714 

.1362 

¥' 

30 

.5888 

.5557 

,5218 

.4885 

.4557 

.4232 

•3904 

.3582 

.3266 

•2957 

.2622 

5- 

■r 

i  40 

.6866 

.6605 

.6350 

•6097 

.5843 

.5589 

.5336 

.5085 

-4837 

.4591 

.4348 

I  50 

.7463 

.7253 

.7046 

.6842 

.6641 

•6433 

.6227 

.6022 

.5819 

.5618 

.5416 

| 

!  75 

.8295 

.8151 

.6011 

.7873 

.7736 

•7594 

.7433 

•7313 

.7168 

.7028 

.6669 

*f„ 

100 

.8707 

.8597 

.8491 

.8385.. 

.8261 

.8172 

.8064 

•7957 

.7830 

.7744 

.7638 

ic 

v. 

I  150 

.9129 

.9054 

.8982 

.8910 

.8840 

.8766 

.8692 

.8619 

.8546 

.8473 

.8401 

« 

??, 

«■ 

I  200 

.9343 

.9887 

.  .9232 

.9178 

.9124 

.9068 

.9012 

.8957 

.8902 

.8847 

.8792 

f 

'  300 

.9559 

.9522 

•9485 

.9448 

.9413 

•9375 

.9337 

.9300 

.9262 

.9825 

.9188 

w 

4oo 

I 

.9669 

.9640 

.9612 

.9585 

.9558 

.9530 

.9501 

.9473 

.9445 

.9417 

.9389 

St 

t 

500 

•9734 

.9712 

.9689 

.9667 

.9646 

.9623 

.9600 

•9577 

•9555 

.9532 

.9510 

^mmmrnw, mt»xw  •; 
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V&Ime  of 


TABLE  4(e)  -  coot. 

1 _ _ 

1  +  (£§)  Fa(2r^2,2n-2p) 


for  a  m  .50 


10  11 


12 


13 


1* 


15 


16 


17 


18 


19 


20 


11 

.061* 

12 

.1268 

.0*i62 

13 

.2016 

.1266 

.0521 

Ik 

.2572 

.1866 

.1171 

.0*&> 

w15 

.3062 

.2392 

.17*7 

.1096 

.0*51 

16 

.3*8* 

.2860 

.2317 

.1630 

.1031 

.0*2* 

17 

.3865 

.3289 

.2681 

.2123 

.1538 

•0973 

.0*00 

18 

.*186 

.3638 

•309* 

.2538 

.1995 

.1*56 

.0921 

.0379 

19 

,**90 

.3976 

•3*55 

.2938 

.2*10 

,109* 

.1382 

.087* 

.0360 

20 

.*757 

.4266 

.3772 

.3268 

.2797 

•2277 

.1803 

.1316 

.0632 

.03*2 

25 

.5791 

.5*00 

-5000 

■*599 

.*207 

■3799 

.3*62 

.2935 

.2635 

-2239 

30 

.6*86 

.61-56 

.5826 

•5*95 

.516* 

•*836 

.*511 

.*186 

.3838 

•351* 

ko 

.7356 

.7109 

.6861 

.6612 

.6363 

.6117 

.5870 

.5622 

•5375 

.5125 

50 

.7882 

.768* 

.7*85 

.7286 

.7085 

.6888 

.6690 

.6*92 

.6293 

.6095 

75 

.8585 

.8*52 

.8320 

.8186 

.8033 

.7921 

.7789 

•7657 

•752* 

•7391 

100 

.8937 

.8838 

.8738 

.6638 

.8538 

.8*39 

.83*0 

.82*0 

.81*1 

.80*1 

150 

•9291 

.922* 

.9158 

.9091 

.902* 

.8958 

.6892 

.8826 

.8759 

.8693 

200 

.9*88 

.9*18 

.9368 

•9318 

.9268 

.9218 

.9168 

.9119 

.9069 

.9019 

300 

.96*5 

.9612 

•9579 

-95*5 

•  9512 

.9*79 

«9**5 

.9*12 

.9379 

•93*6 

*00 

.973* 

.9709 

.968* 

.9659 

•963* 

•  9609 

.958* 

•9559 

-953* 

.9509 

500 

.9787 

.9767 

.97*7 

•9727 

.9707 

.9687 

.9667 

.96*7 

.9627 

.9607 

£ 
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t 


20  item*  are  drawn  at  random  from  a  lot  and  placed  on  life  test. 

The  test  rune  for  100  hours.  Suppose  that  no  failures  occur,  then  we 
can  be: 

(a)  95#  confident  of  the  assertion  that  at  leaet  79,4$  of  the  items 

in  the  lot  survive  100  hours, 

(to)  95 jt  confident  of  the  assertion  that  at  least  86,1%  of  the  items 

in  the  lot  survive  100  hours, 

(c)  90%  confident  of  the  assertion  that  at  least  &9*1£  of  the  items 
in  the  lot  survive  100  hours, 

(d)  75%  confident  of  the  assertion  that  at  least  93. 3$  of  the  items 
in  the  lot  survive  100  hours, 

(e)  «Q%  confident  of  the  assertion  that  at  least  96.6%  of  the  items 
in  the  lot  survive  .100  hours. 

'these  are  non  parametric  statements. 
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3. 

TA2LB  9(«) 


▼aluaa  of  »  ■■-'  ;■  . .  —  l-  for  01  «  .01 

1 4  _ _ _ 


0 

r 

1000 

5,000 

10,000 

50,000 

100,000 

500,000 

1,000,000 

0 

•9954a 

.99917 

•99954 

.99992 

•-99995 

999992 

.999995 

1 

•9933S 

.99867 

•99933 

.99987 

•99993 

.999987 

.999993 

2 

.99162 

.99832 

99916 

.99983 

.99992 

•999983 

.999992 

3 

.98999 

•99799 

.99899 

.99980 

99990 

.999980 

.999990 

4 

.98844 

•99767 

.99684 

•99977 

•99988 

•999977 

.999988 

5. 

.98696 

•99738 

.99669 

.99974 

.99987 

•999974 

.999987 

6 

.98549 

.99708 

.99854 

•99971 

.99985 

•999971  ’ 

.999985 

7 

•96405 

•99679 

.99839 

.99968 

.99984 

•999968 

•999984 

8 

.96270 

.99651 

.99826 

•99965 

.99983 

•999965 

.999983 

9 

.98129 

.99623 

.99811 

.99962 

•99981 

.999962 

.999981 

10 

•97997 

.99596 

•99798 

.99960 

.99980 

•999960 

.999980 

11 

.97863 

•99569 

•99784 

•99957 

.99978 

•999957 

.999978 

12  ' 

•97730 

•99542 

•99771 

.99954 

•99977 

•999954 

•999977 

13 

.97605 

.99517 

-99753 

•99952 

•99976 

•999952 

•999976 

14 

-97465 

.99488 

.99744 

.99949 

.99974 

.999949 

•999974 

15 

•97334 

.99462 

•99731 

•99946 

•99973 

•999946 

•999973 

16 

•97209 

.99436 

.99718 

.99944 

.99972 

•999944 

•999972 

17 

.97085 

.  .99411 

•99705 

.99941 

•99971 

.999941 

•999971 

18 

.96961 

.99386 

.99693 

.99938 

.99969 

•999938 

.999969 

19 

.96841 

•99362 

.93680 

.99936 

•99968 

•999936 

•999968 

20 

.96717 

.99336 

•99668 

•99933 

•99967 

•999933 

•999967 

30 

.95490 

.99087 

.99543 

.99908 

•99954 

-999908 

•999954 

40 

.94318 

.98846 

.99423 

.99884 

.99942 

.999884 

•999942 

50 

-93149 

.96618 

.99308 

.99860 

•99931 

.999860 

•999931 

60 

.91962 

.98364 

.99180 

.99836 

•99918 

.999836 

•999918 

70 

•90807 

.98126  . 

.99061 

.99812 

•99906 

-999812 

.999906 

80 

.89695 

.97897 

.98946 

.99789 

.99894 

.999789 

.999894 

90 

.88621 

.97676 

.98835 

.99766 

.99883 

•999766 

.999883 

100 

.87368 

.97488 

.96726 

.99745 

.99874 

.999745 

.999874 

200 

.76750 

.95267 

.97628 

.99525 

•99762 

•999525 

•999762 

500 

.46461 

.89242 

.94616 

.98923 

•99461 

.998923 

•999461 

j. 
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MBLi  9(b) 


Values  of 


1 

1  *  {~i)r  (2r^.2n.2r) 

n*3r  U 


for  o  -  .05 


- T 

r 

1000 

5000 

10.000 

50,000 

— 

100,000 

500,000 

100,000 

0 

.3? 

.99703. 

^999*0 

•99970 

•9999* 

.99997 

qqqqqIi 

*77777* 

.999997 

1 

v^527 

•99905 

-99953 

.99991 

•99995 

.999991 

qooooh 

2 

•99371 

•9967* 

•99937 

,99987 

•9999* 

.999987 

.99999* 

3 

.99226 

.990*5 

.99922 

•9998* 

•99992 

.99998* 

.999992 

It 

•99007 

•99017 

,99906 

.99902 

99991 

■>999982 

.999991 

5 

•96953 

•99790 

•  99095 

•99979 

.99989 

.999979 

•999989 

6 

.96620 

.99763 

,99881 

•99976 

.99988 

•999976 

.999988 

7 

•90692 

•99737 

,99868 

•99974 

.99987 

•  99997* 

.999987 

6 

•90565 

•99711 

.99856 

-99971 

99986 

.999971 

.999986 

9 

,96*36 

-99685 

,998*3 

.99969 

•9998* 

•999969 

.99996* 

10 

.96312 

.99661 

.99830 

•99966 

•99963 

.999966 

.999983 

11 

■96163 

•99635 

.99817 

-99963 

•99982 

.999963 

•999982 

12 

.98058 

•  99609 

.9900* 

.99961 

<99900 

,999961 

•999900 

13 

•97937 

•99505 

•9^792 

.99958 

•99979 

.999958 

•999979 

1* 

.97020 

•  99561 

•99760 

.99956 

,99978 

•999956 

•999978 

15 

.97697 

•99536 

-99768 

•99954 

•99977 

,99995* 

•999977 

16 

97570 

.99512 

•99756 

99951 

,99976 

*999951 

•999976 

17 

97*59 

•99*06 

.9974* 

•999*9 

•9997* 

•9999*9 

,999974 

16 

•973*1 

•  99*6* 

•99732 

•999*6 

-99973 

.9999*6 

.999973 

19 

•97225 

0 99**1 

,99720 

•999*4 

•99972 

.9999** 

,999972 

20 

•97106 

.99*17 

,99706 

•999*2 

•99971 

.9999*2 

-999971 

30 

•95951 

.99183 

•99591 

•99916 

•99959 

.999918 

•999951 

*0 

.9*626 

.90955 

•9r*77 

.99895 

.999*8 

999895 

•9999*8 

50 

-93711 

.98738 

.  99369 

.99874 

.99936 

.99987* 

•999936 

60 

<96503. 

,98*98 

,99248 

•990*9 

■  .99925  : 

.9998*9 

•999925 

70 

.91*76 

•98273 

•99135 

.99827 

99913 

.999827 

•999913 

do 

•90397 

.98052 

,9902* 

.99805 

.99902 

.999605 

-999902 

90 

•69337 

.97852 

.90925 

99785 

.99692 

999785 

100 

.88297 

•976*1 

,98819 

-9976* 

.99882 

*990764 

•999082 

200 

•77703 

.95*  f* 

•97733 

•995*6 

.99773 

.9995*6 

•999773 

500 

.*7456 

.99*39 

.9*717 

•989*3 

.99*71 

.9989*3 

* 

•999*71 

Value*  of 


for  a  »  .10 


3.  ‘>W 


TABLE  Me) 

1  _ 

l*  <~>  ro( 2r*2,2u-2r) 


f 


a 

r 

1000 

5000 

10,000 

50,000 

100,000 

500,000 

1,000,000 

0 

•99771 

.99954 

.99977 

•99995 

.99998 

.999995 

.99999/ 

1 

•99612 

•99922 

•99961 

.99992 

.99996 

.999992 

.999/6 

2 

•  99469 

.99094 

.99947 

.99989 

.99993 

.999989 

.956995 

3 

•99333 

.99866 

•99933 

.99987 

.99993 

.999987 

,*<999993 

4 

.99202 

.99840 

•99920 

■  99964 

.99992 

.999984/ 

•999992 

*; 

•99069 

•99013 

.99907 

.99981 

.99991 

&99^1 

•999991 

6 

•90945 

•99706 

•  99894 

•99978 

.99989 

/.  999979 

.999989 

7 

•98822 

,99764 

.99882 

•99976 

.99988 

■999976 

•999988 

8 

•96676 

.99734 

.99667 

.99973 

.99987 

•999973 

.999987 

9 

.9058? 

•99716 

.99858 

•99971 

.99986 

•999972 

•999986 

10 

•98459 

.99691 

•  99846. 

.99969 

.99985 

.999969 

.999985 

11 

•90329 

.99666 

.99834 

-99968 

.99983 

•999968 

•999983 

I 

12 

,90205 

•99642 

.99823 

•99966 

.99982 

•999966 

•999982 

13 

.98087 

•99617 

.99811 

•99965 

.99981 

•999965 

.999981 

14 

-97970 

•99593 

•99799 

•99964 

.99980 

•999964 

•999900 

15 

•97051 

•99570 

.99787 

•99962 

.99979 

•  9999:2 

•999979 

16 

.97736 

•99547 

•99776 

•99961 

.99978 

.999961 

.999978 

17 

.97621 

.99523 

•99764 

.99960 

.99977 

•999960 

.999977 

18 

•97509 

•99501 

•99752 

.99958 

.99975 

■999958 

.999975 

19 

•97401 

.99479 

-99740 

•99957  , 

.99974 

•999957 

.999974 

20 

•97296 

.99450 

•99729 

■99956  ' 

rf973 

•999946 

•999973 

30 

.96154 

•99231 

•99615 

•99923' 

•999923 

.999962 

40 

•95131 

.99002 

.99513 

•99900 

.99951  - 

•999900 

.999951 

50 

•94010 

.98798 

<99399 

.99880 

.99940 

.999880 

•999940 

60 

•92937 

.90563 

•99294 

•99056 

.99930 

•999856 

•999930 

70 

.91876 

•90342 

.99180 

.99834 

.99919 

.999834 

•  999*319 

00 

.90762 

.90136 

.99076 

.99813 

.99908 

.999613 

•999908 

90 

.89693 

.97929 

.98969 

.99793 

.99897 

•999793 

.9^9897 

100 

.88675 

•97721 

.90660 

•99772 

.99886 

.999772 

•999006 

200 

.78227 

.95603 

•97799 

•99559 

.99700 

•999559 

•999700 

500 

.40201 

.89591 

•94792 

.98958 

.99479 

.998958 

.999479 

1000 


1000 

5000 

10,000 

0 

•99972 

.99906 

1 

•99730 

•99946 

.99973 

2 

•99607 

•99921 

•99961 

3 

•99489 

.99696 

.99949 

4 

•99376 

.99875 

•99938 

5 

-99257 

.99651 

•99926 

'  6 

•99143 

•99028 

-99914 

7 

.99030 

•99006 

■99903 

8 

.90920 

•99784 

.99892 

9 

.96812 

•99762 

P9S81 

10 

.98698 

-99739 

.99070 

11 

.90584 

•99716 

*99850 

12 

.98472 

•99693 

•99847 

13 

.  98.362 

•99671 

•99036 

14 

.98254 

,.99648 

.99825 

15 

•90146 

.99626 

.99015 

16 

.96030 

.99605 

.99603 

17 

•97926 

•99584 

•99793 

18 

-97022 

•99563 

19 

•97718 

-99542 

•  '99772 

20 

••97614 

•99522 

,99761 

30 

.96459 

,99308 

•99646 

40 

■95439 

•99096 

.99544 

50 

.94448 

.96887 

.99443 

60 

.93445 

.98684 

•99344 

70 

.92443 

.90481 

•99243 

80 

.91441 

.96282 

•99143 

90 

.90440 

.98063 

.99044 

100 

.89439 

•97805 

.90942 

200 

.79437 

.95005 

.97942 

500 

.49434 

.89664 

.94942 

•999997 


•99990 
•99997 
*99996 
•99995 
•99994 
•99993 
•99991 
.99990 
<99909 
•99968 
•  99987 
•99906 
•99905 
.99904 
•99903 
•99902 
•99900 
•99979 
•99978 
•99977 
•99976 
•99965 
•99954 
.99944 
.99934 
•99924 
.99914 
.99904 
•99094 
•99794 
•99494 


■999972 
•99996 9 


•999952 

•999931 

•999910 


•999995 


•999909 

•999908 

•999907 

•999906 

*999905 

•999904 

*999903 

•999902 

■999900 

•999979 


TABU  5(e) 


Values  of  -- —■  r  -  "  »  for  a  «  .50 

1  ♦  (—)  r(2r*2,2n-2r: 
n-r  0 


^3 

1000 

5000 

10, >300 

50,000 

100,000 

500,000 

1,000,000 

0 

•99931 

•99986 

•99993 

■99999 

99999 

999999 

.999999 

1 

•99832 

•99966 

.99983 

•99997 

99998 

.999997 

999998 

2 

•99733 

•99947 

•99973 

•99995 

99997 

-999995 

-999997 

3 

•99633 

.99927 

99963 

•99993 

•  99996 

•999993 

.999996 

4 

•99533  " 

■.99907 

•99953 

•99991 

99995 

999991 

999995 

5 

.99435 

.99887 

.99943 

.99989 

.99994 

•999989 

*999994 

6 

•99334 

•99867 

■99933 

.99987 

■99993 

999987 

•999993 

7 

.99234 

.99847 

.99923 

.99985 

*99992 

999985 

999992 

8 

•99134 

.99827 

.99913 

■99963 

•99991 

.999983 

•999991 

9 

•99033 

.99807 

•99903 

.99961 

99990 

•999981 

•999990 

10 

•98934 

■99787 

.99893 

•99979 

99989 

■999979 

.999989 

u 

.98832 

•99767 

•99883 

.99977 

■99988 

•999977 

.999988 

12 

•98730 

99747 

.99873 

>99975 

.99987 

•999975 

.999987 

13 

.98630 

.99728 

.99863 

.99973 

.99986 

•999973 

•999986 

14 

.98531 

.99708 

.99853 

•99971 

•99985 

• 999971 

.999985 

15 

•98432 

.99687 

•99843 

.99969 

99984 

.999969 

•999984 

16 

.98333 

.99666 

.99833 

•99967 

99983 

.99.99=7 

.999983 

17 

.98234 

•99646 

.99823 

99965 

.99982 

.999965 

.999982 

18 

.98135 

.99627 

.99814 

.99963 

.99981 

*999963 

.999981 

19 

.98036 

.99607 

-99804 

•99961 

.99980 

•999961 

999980 

20 

.97936 

.99587 

.99793 

99959 

-99979 

♦999959 

•999979 

30 

.96937 

.99386 

•99693 

•99939 

■99969 

-999939 

•999969 

40 

.95936 

•99X86 

-.99593 

•99919 

*99959 

999919 

•999959 

50 

.94936 

.98987 

.99494 

.99899 

.99949 

.999899 

999949 

60 

•93935 

.98768 

99394 

.99879 

*99939 

.999879 

*999939 

70 

.92934 

.98588 

.99293 

.99859 

•  99929 

.999859 

*999929 

80 

•91935 

.98387 

.99193 

.99839 

99919 

999839 

•999919 

90 

.90936 

.98186 

.99094 

■99819 

.99909 

.999819 

•999909 

100 

.89935 

■97986 

.98993 

•99799 

99899 

•999799 

*999699 

boo 

•79942 

.95987 

-97993 

•99599 

•99799 

•999599 

-999799 

J500 

.49956 

.89986 

.94993 

*98999 

*99499 

.996999 

.999499 

ftaMr.'  cal  Example  far  SHu&le*  y 


1000  Items  are  drain  at  random  from  a  lot  sM  plated. -on  life  test-  'Jhe 
test  runs  for  100  hour-  Suppose  that  20  faiiur*i?  cx-eur,.  then  we  can  be : 


(a)  993>  co.ifident  of  the  assertion  that  at  least  96  72^  of  the 
items  in  the  let  survive  .100  hours: 

(b)  95^  confide  V-  cl’  i^swrttoa  that  at  least  97 «■  li  $  as'  tha 

items  in  th  t  J*i  ACC*  K  >urr ; 

(c)  90 ft  confid:  «i  of  the  -  that  at  least  97. 30^  of  the 

Items  in  tbt.  Ice  if.  rvi*:-  100  hour;; 

(d)  75jt  confident  of  the  &»s»2ition  that  at  least  97  5.1  of  the 

items  In  the  lot-  survive  100  hours;  * 

(•)  50^  confident  of  the  assertion  that  at  leant  Q?  94jt  of  the 


Items  in  the  lot  i-urvive  100  hour... 


